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Abstract 

Ghost  imaging  is  a  transverse  imaging  technique  that  relies  on  the  correlation  between  a 
pair  of  light  fields,  one  that  has  interacted  with  the  object  to  be  imaged  and  one  that  has 
not.  Most  ghost  imaging  experiments  have  been  performed  in  transmission,  and  virtually  all 
ghost  imaging  theory  has  addressed  the  transmissive  case.  Yet  stand-off  sensing  applications 
require  that  the  object  be  imaged  in  reflection.  We  use  Gaussian-state  analysis  to  develop 
expressions  for  the  spatial  resolution,  image  contrast,  and  signal-to-noise  ratio  for  ghost 
imaging  performed  by  measuring  a  fraction  of  the  light  that  reflects  off  a  rough-surfaced 
object  that  creates  target  returns  with  fully  developed  speckle.  This  is  done  for  a  pseu- 
dothermal  source  with  phase-insensitive  classical  correlation  between  the  two  fields,  and  for 
a  quantum  source  with  non-classical  phase-sensitive  correlation  between  the  fields.  In  the 
low  flux  limit,  this  quantum  source  becomes  the  biphoton  state.  We  compare  our  results  to 
the  corresponding  behavior  seen  in  transmissive  ghost  imaging,  and  we  develop  performance 
results  for  the  reflective  form  of  computational  ghost  imaging.  We  also  provide  a  preliminary 
stand-off  sensing  performance  comparison  between  reflective  ghost  imaging  and  a  conven¬ 
tional  direct-detection  laser  radar  system.  We  also  consider  the  resolution  degradation  on 
each  system  when  the  fields  propagate  through  turbulence.  Finally,  we  investigate  ways  of 
increasing  the  signal-to-noise  ratio  of  reflective  ghost  imaging  through  use  of  multiple  bucket 
detectors,  multiple- wavelength  sources,  and  compressive  sensing. 
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Chapter  1 


Introduction 


In  this  thesis  we  explore  the  viability  of  using  ghost  imaging  for  stand-off  sensing.  To 
this  end  wc  characterize  the  imago  quality  of  such  systems,  and  compare  the  results  to  the 
image  quality  of  an  equivalent  laser  radar  system.  Wc  connect  our  work  back  to  previous 
theory  developed  for  transmissive  ghost  imaging,  as  well  as  to  experiments  done  in  reflection. 
Finally,  we  explore  ways  to  improve  the  quality  of  reflective  ghost  imaging.  The  most 
promising  of  these  approaches  is  the  application  of  compressive  sensing,  which  leverages  the 
structure  of  natural  images  to  reduce  image  capture  time  while  improving  image  quality. 

Ghost  imaging  exploits  the  correlation  between  two  light  fields  to  create  an  image  that 
neither  field  alone  could  provide.  We  denote  these  two  fields  the  ‘'signal”  and  “reference" : 
the  signal  field  interacts  with  the  target,  after  which  a  single-pixel  “bucket"  detector  makes 
a  power  measurement  of  the  field;  simultaneously,  the  reference  field’s  transverse  power  dis¬ 
tribution  is  measured  with  a  high  spatial-resolution  detector,  which  is  usually  a  scanning 
pinhole  or  a  CCD  array.  The  signal  and  reference  fields  have  some  cross  correlation,  de¬ 
termined  by  the  choice  of  source,  so  measuring  the  reference  field  provides  some  knowledge 
of  the  field  illuminating  the  target.  Since  the  power  measurement  on  the  signal  arm  is  a 
function  of  the  target  and  the  field  illuminating  it.  we  are  able  to  reconstruct  the  target  by 
correlating  the  power  measurement  with  the  output  of  the  high  spatial-resolution  detector. 
This  imaging  technique  has  become  known  as  ghost  imaging  because  the  image  information 
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is  contained  in  the  correlation  between  the  two  measurements:  neither  measurement  alone 
contains  sufficient  information  to  reconstruct  the  target. 

The  first  ghost  imaging  experiment  used  biphoton  pairs  for  the  signal  and  reference 
fields,  which  were  generated  by  spontaneous  parametric  downconversion  (SPDC)  and  post¬ 
selection  [1].  Because  biphotons  are  entangled  states—  with  a  phase-sensitive  cross  correlation 
between  the  signal  and  reference  fields  stronger  than  allowed  by  classical  physics  Tor  which 
the  quantum  theory  of  photodetcction  is  needed  to  calculate  the  measurement  statistics,  it 
was  initially  thought  that  ghost  imaging  wa|  a  non-local  quantum  phenomenon.  However, 
the  non-classical  nature  of  ghost  imaging  was  called  into  question  when  experiments  were 
later  performed  using  a  classical  pseudothermal  light  source  [2.  3].  In  this  setup  laser  light 
is  passed  through  a  ground-glass  diffuser,  after  which  it  is  divided  by  a  50  -50  beam  splitter 
into  identical  signal  and  reference  fields  with  a  phase-insensitive  cross  correlation. 

Gaussian-state  analysis  has  enabled  a  unified  treatment  of  biphoton  and  pseudothermal 
ghost  imaging  that  shows  the  image  formation  process  is  one  of  classical  coherence  propaga¬ 
tion,  with  high  contrast  in  DC-coupled  biphoton  ghost  images  being  the  principal  ghost-image 
signature  of  that  non-classical  source  [4],  In  particular,  identical  statistics  for  pseudothermal 
ghost  imaging  result  from  the  use  of  quantum  photodetectiou  theory,  in  which  both  the  light 
beams  and  photodetectors  arc  treated  quantum-mechanically;  and  scmiclassical  photodetec¬ 
tion  theory,  in  which  the  light  beams  are  treated  classically  but  photodotcctors  inject  shot 
noise  on  top  of  any  fluctuations  in  the  illumination.  A  thorough  review  of  these  considera¬ 
tions  can  be  found  in  [5].  For  our  purposes,  it  suffices  to  note  that  we  can-  and  will  -use 
scmiclassical  photodetection  theory  when  dealing  with  classical  sources. 

To  date,  Gaussian-state  analysis  of  transmissive  pseudothermal  ghost  imaging  has  pro¬ 
vided  expressions  for  its  spatial  resolution,  image  contrast  and  signal-to-noise  ratio  (SNR) 
behaviors;  i.e..  a  complete  characterization  of  its  performance  [4,  6].  However,  one  of  the 
more  interesting  potential  applications  of  ghost  imaging  is  stand-off  sensing,  in  which  the 
bucket  detector  observes  the  target  in  reflection,  not  transmission.  Preliminary  table-top  ex¬ 
periments  have  demonstrated  the  feasibility  of  this  approach  [7,  8] ,  but  there  has  been  little 
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exploration  of  the  statistical  characteristics  of  these  images.  Developing  that  theory  within 
the  Gaussian-state  framework,  which  is  a  focus  of  this  thesis,  must  confront  an  additional 
complication  not  seen  in  previous  work:  viz.,  the  speckle  induced  by  reflection  from  rough- 
surfaced  objects.  In  particular,  we  shall  report  expressions  for  the  spatial  resolution,  image 
contrast,  and  SNR  of  reflective  ghost  imaging  of  speckle  targets  for  classical  and  non-classical 
sources  in  Chapters  2  and  3,  respectively.  These  will  be  compared  with  the  corresponding 
results  for  the  transmissive  case,  in  which  there  is  no  target-induced  speckle.  In  Chapter  4 
we  will  use  this  framework  to  obtain  performance  results  for  the  reflective  form  of  computa¬ 
tional  ghost  imaging  [9,  10],  in  which  the  reference  beam  is  removed.  Then,  in  Chapter  5, 
we  will  provide  a  preliminary  stand-off  sensing  performance  comparison  between  reflective 
ghost  imaging  and  a  conventional  direct-detection  laser  radar. 

The  effects  of  turbulence  on  transmissive  ghost  imaging  have  been  explored  experimen¬ 
tally  [8]  and  theoretically  [11].  As  most  practical  applications  of  reflective  ghost  imaging  will 
require  propagation  through  atmospheric  turbulence,  in  Chapter  6  we  extend  our  analysis 
to  include  turbulence  on  all  paths.  The  effects  of  turbulence  on  image  resolution  will  bo 
explored  for  classical,  non-classical,  and  computational  reflective  ghost  imaging,  as  well  as 
for  a  laser  radar  system.  In  Chapter  7  we  explore  ways  to  improve  the  image  quality  of 
reflective  ghost  imaging.  The  first  of  these  methods  is  to  employ  multiple  c.o-planar  bucket 
detectors  on  the  signal  arm;  since  the  target-induced  speckle  decorrelates  with  transverse 
separation  in  the  detector  plane,  each  bucket  detector  observes  different  speckle  behavior, 
allowing  us  to  average  out  its  deleterious  effects  on  the  SNR.  The  second  method  is  to  use 
multiple  sources  at  different  wavelengths.  This  again  seeks  to  average  out  the  target-speckle 
effects:  if  the  wavelengths  arc  far  enough  apart,  the  returns  associated  with  each  source  have 
uncorrelated  speckle  statistics.  Finally,  we  extend  the  compressive  sensing  methods  used  in 
transmissive  imaging  [12]  to  reflective  imaging.  Compressive  sensing  allows  for  reduced  im¬ 
age  capture  time  and  higher  quality  reconstructions  when  the  target  is  sparsely  represented 
in  some  basis  [13.  14],  Chapter  8  concludes  the  thesis  with  a  brief  summary  of  results  and 
suggestions  for  further  research. 
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Chapter  2 


Ghost  Imaging  with  Pseudothermal 
Light 

2.1  Theoretical  Setup 

In  stand-off  sensing,  targets  will  most  likely  be  separated  from  the  source  by  distances  on  the 
order  of  kilometers.  This  indicates  that  high-powered  classical  sources  might  be  necessary. 
While  we  also  consider  low-flux  quantum  sources  in  Chapter  3,  we  start  our  analysis  with 
the  system  shown  in  Fig.  2-1,  which  utilizes  pseudothermal  light  to  perform  reflective  ghost 
imaging.  A  continuous-wave  (cw)  laser  beam  is  passed  through  a  rotating  ground-glass 
diffuser  followed  by  a  50  50  beam  splitter  to  produce  identical,  spatially-incoherent  signal  and 
reference  beams  whose  temporal  bandwidths  are  much  lower  than  those  of  the  single-pixel 
(bucket)  and  high  spatial-resolution  (CCD  array)  detectors.  The  signal  beam  illuminates 
a  rough-surfaced  planar  target  at  distance  L  from  the  beam  splitter,  and  some  of  the  light 
reflected  from  that  target  is  collected,  after  L-m  propagation,  by  the  bucket  detector.  The 
reference  beam  directly  illuminates  the  CCD  array  which,  for  theoretical  convenience,  we 
have  placed  L-m  away  from  the  beam  splitter.1  The  photocurrents  from  the  bucket  detector 

1This  assumption  implies  that  we  will  form  a  1:1  ghost  image.  In  an  actual  implementation  of  reflective 
ghost  imaging  the  CCD  array  would  be  in  the  focal  plane  of  a  lens  located  near  the  transmitter  and  we 
would  obtain  a  minified  ghost  image. 
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Figure  2-1:  Setup  for  pseudothermal  reflective  ghost  imaging. 


and  each  pixel  on  the  CCD  are  sent  to  a  correlator,  whose  output  for  the  CCD  pixel  located 
at  transverse  coordinate  pi  is 


i  rTi'2 

C(pi)  =  —  /  dtii{t)i2{t).  (2.1) 

1 1  J-Ti/2 


where  Tj  is  the  averaging  time  and  we  have  suppressed  an  L/c  time  delay  in  ii(t)  that  is 
needed  to  account  for  the  delay  incurred  by  the  light,  returning  from  the  target. 


The  configuration  and  notation  we  are  using  parallels  the  semiclassical  treatment  of 
transmissive  ghost  imaging  in  [9]  with  the  principal  distinction  being  that  in  the  transmissive 
case  the  bucket  detector  would  be  behind  a  transmission-mask  target  L-m  from  the  signal 
source,  whereas  here  that  target  is  viewed  in  reflection.  We  have  also  switched  the  labels 
for  the  signal  and  reference  fields;  in  [9]  the  reference  field  illuminates  the  target,  while  the 
signal  field  goes  to  the  high  spatial-resolution  detector. 


►  Correlator  A 


hit) 


All  the  fields  shown  in  Fig.  2-1  are  complex  envelopes  about  center  frequency  cuq  of 


linearly-polarized  light  fields  normalized  to  have  v/photons/m2s  units  as  functions  of  their 
transverse  coordinate  vectors  and  time.  As  a  result,  under  the  assumption  of  shot-noise 
limited  detectors  with  quantum  efficiency  77,  the  photocurrents  from  the  bucket  and  the 
CCD  can  be  written  as 

im(t)  -  J dr  [ qr]Pm{T )  +  A im(r)]  h{t  -  r),  for  m  =  1.2.  (2.2) 

Here:  PJi)  =  JA  dp  \Em(p,  t)\2  is  the  photon  flux  on  the  circular  active  region  Am  of 
detector  rn\  q  is  the  electron  charge;  A im(t)  is  the  shot  noise  from  detector  m;  and  h(t) 
is  the  detector’s  baseband  impulse  response,  normalized  to  satisfy  Jdth(t)  =  1.  Phys¬ 
ically,  qi]Pm(t)  is  the  conditional  mean  of  given  the  illumination,  so  that  A im(t) 

is  the  photocurrent  fluctuation  conditioned  on  knowledge  of  the  illumination.  Note  that 
given  the  photon-flux  waveforms  Pi(t)  and  Po(t),  the  shot-noise  currents  Aij (t)  and  A i2{t) 
are  statistically  independent,  zero-mean,  random  processes  whose  correlation  functions  are 
(Afm(n)Afm(t2))  =  q2r]Pm{ti)6(t1  -  to). 


The  fields  that  determine  the  preceding  photon  fluxes  are  found  from  diffraction  theory, 
applied  on  three  separate  paths:  the  reference  path  (R).  from  the  source  to  the  high  spatial- 
resolution  detector;  the  signal  path  (S),  from  the  source  to  the  target:  and  the  target-return 
path  (T),  from  the  target  to  the  bucket  detector.  I11  the  ensuing  work,  we  will  denote  fields 
that  have  propagated  L- m  on  path  m  with  a  prime;  i.e.,  if  wc  start  with  some  field  E  m(p'P), 
after  L-m  propagation  it  becomes  E'm(p,t )  as 

f  bn»Po{L+\p-p'[2pL) 

E'm(p ,  t)  =  j  dp'  Em(p',  i)- - — -  m=R,  S.  T.  (2.3) 

where  k0  =  Uq/c  is  the  wave  number,  and  we  have  suppressed  time  delays.  The  fields 
are  defined  as  follows:  ER(p.t )  is  the  reference  field  at  the  source,  and  E'R{p,t)  is  the 
field  illuminating  the  CCD;  Es{p.t)  is  the  signal  field  at  the  source,  and  E's{p,t)  is  the 
field  illuminating  the  target;  ET{pp)  —  Ers(p,t)T(p)  is  the  target-return  field,  which  is 
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the  propagated  signal  field  immediately  after  reflecting  off  a  target2  with  field-reflection 
coefficient  T(p),  and  E'T(p,t)  is  the  field  illuminating  the  bucket  detector.  Since  we  are 
using  the  semielassical  theory  of  light,  each  photodetector  is  taken  to  directly  measure  the 
classical  field  impinging  it,  so  Ei(p.  t)  —  E'R{p,  t),  and  E&p.t)  =  E'r(p.  /.). 

Whereas  in  transmissive  ghost  imaging  it  is  ordinarily  the  case  that  the  target’s  field- 
transmission  coefficient  is  assumed  to  be  deterministic,  the  targets  of  interest  for  reflective 
ghost  imaging  will  have  microscopic  surface  variations  from  a  nominal,  smooth  surface 
profile  whose  standard  deviations  can  greatly  exceed  the  illumination  wavelength  and  whose 
transverse  correlation  scale  can  be  sub-wavelength.  When  such  a  surface  is  illuminated  by 
laser  light  it  gives  rise  to  laser  speckle  in  the  target  return,  and  a  reasonable  statistical  model 
for  that  behavior  is  to  take  the  target-  reflectivity  to  be 

T(p)  =  ^T{p)e2,knAz(p),  (2.4) 

where  the  height  Variations  A z(p)  are  modeled  as  a  zero-mean,  real- valued,  Gaussian  random 
process  with  a  transverse  correlation  on  the  order  of  a  wavelength.  Thus  T(p)  is  a  zero-mean, 
complex- valued  Gaussian  random  process  that  is  completely  characterized  by  the  correlation 
function  [15] 


(T*(Pl)T(P2))  =  A 20T(Pl)6(Pl  -  p2),  (2.5) 

where  A0  is  the  center  wavelength  of  the  illumination  and  T(p)  is  the  target's  intensity- 
reflection  coefficient,  which  is  nonrandom  and  the  quantity  that  we  are  seeking  to  image 
with  the  Fig.  2-1  setup. 

In  order  to  proceed  further,  we  need  to  specify  the  source-field  characteristics.  Following 
the  Gaussian-state  analysis  of  Erkmen  and  Shapiro  [4,  6,  9],  we  shall  assume  that  Es(p,  t) 
and  Er(p,  t) — the  identical  outputs  from  the  50-50  beam  splitter  in  Fig.  2-1 — are  zero-mean, 
complex-valued  Gaussian  random  fields  that  are  completely  characterized  by  their  common 

2We  have  assumed  a  stationary  target,  so  that  its  field-reflection  coefficient  is  constant  in  time. 
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cross- spectrally  pure  phase-insensitive  correlation  function  per  the  Gaussian-Schell  model  as 


{El{PlM)Em{p2.t2))  =  ^e-^d2+|P2|2}/ag-|P1-p2|V2p^-(^nF/2T0^  tm  €  (2.6) 


7ra, 


with  photon-flux  P,  e  2  intensity  radius  a0,  coherence  length  p0,  and  correlation  time  T0. 


The  preceding  setup  fully  specifies  all  that  is  needed  to  compute  the  spatial  resolution, 
image  contrast,  and  SNR  of  the  pseudothermal  reflective  ghost  image.  Before  doing  so, 
however,  let  us  introduce  one  final  condition.  We  shall  assume  that  the  CCD  pixel  active 
region,  .Ax,  is  sufficiently  small  that  \Ei(p,  t)\2  is  essentially  constant  over  each  pixel,  allowing 
us  to  use  Pi(t)  —  Ai\Ei(p.t) |2,  where  A1  is  the  area  of  A\.  This  condition  ensures  that  the 
spatial  resolution  we  obtain  is  limited  by  the  field  statistics,  not  by  the  CCD's  pixel  size. 


2.2  Spatial  Resolution  and  Image  Contrast 


The  spatial  resolution  and  image  contrast  of  the  pseudothermal  reflective  ghost  image  are 
properties  of  the  ensemble-averaged  photocurrent  cross  correlation,  i.e.,  Eq  (2.1)  averaged 
over  the  shot  noise  and  the  fluctuations  in  the  fields  that  illuminate  the  two  detectors.  Those 
field  fluctuations  arise  from  the  randomness  imposed  by  the  ground-glass  diffuser  and,  for  the 
field  illuminating  the  bucket  detector,  the  target’s  surface  roughness,  ft  is  easy  to  see  that 
the  fields  illuminating  the  detectors  are  zero-mean  random  processes  that  are  statistically 
stationary  in  time,  and  thus  from  Eq  (2.2)  the  ensemble- averaged  ghost  image  satisfies 

(npi)>  =  (d(DC(/)> 

=  J drx  J dr2  h(t  -  n)h(t  -  r2)  J  dp  (\ElR(p1.T1)\2\E!r(p' ,  r2)|2),  (2.7) 
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Back  propagating  with  Eq  (2.3)  to  the  field  illuminating  the  target,  the  intensity  correlation 
becomes 


(\E'R(pi.T^\2\E'r{p',T2)\2)  =  J  dp-2  Jdp3eiko{lp’  ^  |p'  P^I2L 

x  (Er{pi.,  T1)E's{p2,T2)E'R(p1,  Ti)Els{p3:  r2))(T*(p2)T(p3)),  (2.8) 


where  we  have  exploited  the  statistical  independence  of  E'm(p,  t)  and  T(p),  i.e..  the  fluctua¬ 
tions  due  to  propagation  through  the  ground  glass  and  those  induced  by  the  target's  surface 
roughness.  Equation  (2.3)  shows  that  E'm(p,t),  for  m  =  R.S,  is  a  linear  transformation 
of  the  zero-mean  Gaussian  random  process  Em(p.t);  hence  it  too  is  a  zero-mean,  Gaussian 
random  process.  Thus,  as  in  the  Gaussian-statc  analysis  of  transmissive  ghost  imaging  [4], 
we  employ  Gaussian  moment- factoring  to  write  the  fourth-order  field  moment  in  terms  of 
second-order  moments.  Making  use  of  Eq  (2.5)  to  evaluate  the  surface  moment,  and  for  A2 
being  the  area  of  A2,  Eq  (2.8)  becomes 


(|£'Upi-r1)|2|£Up',r2)|2) 

=  77  jiPiTM  [<|b«(pi.t1)P)(|£'Up2.%)|2>  +  |( EHPl,n)E’s(P2,T,))f}  (2.9) 

which  clearly  indicates  the  role  of  photon-flux  correlation  in  ghost  image  formation. 

Pseudothermal  ghost  imaging  is  performed  with  spatially  incoherent  light,  i.e.,  p0  -C  a0 
holds  in  the  source's  Gaussian-Schell  model  spatial  correlation.  Moreover,  stand-off  sensing 
pseudothermal  ghost  imaging  will  be  performed  in  the  far  field,  for  which  k0a0p0/2L  <C  1 
prevails.  Following  Erkmen  and  Shapiro  [5],  we  can  propagate  the  Gaussian-Schell  correlation 
function  in  Eq  (2.6)  into  the  far  field  and  obtain 


(E'e*(pi:  1 1  )/•-'(„( Pi>.  h)) 


2  P 


7T  a, 


IPi  l2)p-(IPl  r2+te2i2)/a|-|pl-p2i2/2P'i(O.*-fl2-il)2/2T02 


(2.10) 


for  i.rn  G  {R.S},  with  a i  =  2L/k0po  and  p^  =  2L//co«o  being  the  new  intensity  and 
coherence  radii,  respectively. 
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To  reduce  Eq  (2.7)  to  a  form  in  which  we  can  easily  assess  the  reflective  ghost  image’s 
spatial  resolution  and  image  contrast,  we  need  to  define  detector’s  impulse  response  h(t). 
For  now,  let  us  take  this  to  be  a  DC-coupled  Gaussian  of  bandwidth  i.e., 


h{t)  = 


V-b  nity  8 

V&r 


(2.11) 


Because  we  have  assumed  that  the  detector  bandwidth  is  sufficient  to  follow  the  temporal 
behavior  of  the  light  emerging  from  the  ground  glass,  we  have  that  Q.bTq  >  1.  This  lets  us 
simplify  the  evaluation  of  the  temporal  integrals  and  obtain 


(C'(pi)) 


q2rf  A\Ao  (  2 P 


L2 


Tver, 


dp2  T(p2)e”2(|pl|2+l^|2)/a'  [1  +  (2.12) 


When  the  intensity  radius  aL  is  much  larger  than  the  target’s  transverse  extent,  so  that  the 
entire  target  is  uniformly  illuminated  on  average,  we  get  our  final  form  for  the  ensemble- 
averaged  photocurrent  cross  correlation, 


(C(pi))  = 


q2rj'2  A  i  Ao  f  2 P 


Lr 


7r  a. 


dp->  Tipi)  \  +  e 


-\pi-p->\2  i  p\ 


(2.13) 


Equation  (2.13)  shows  that  the  eiijpnble-  aver  age  photocurrent  cross  correlation  consists 
of  a  featureless  background  term, 


inrA\A2 
L 2 


J  dp2T{p2 ), 


(2.14) 


plus  the  image-bearing  term, 


Ci(pi 


2  2 
q  r) 


AiA2  (  2P\ 


L 2 


/ dp2T(p2y 
nciiJ  J 


-\P\-P'a\2/p2,, 


(2.15) 


Resolution  The  image-bearing  term  contains  the  target’s  intensity-reflection  coefficient 
T(P2)  convolved  with  a  Gaussian  point-spread  function  (PSF)  that  limits  the  spatial  resolu- 

3The  frequency  response  associated  with  this  impulse  response  is  H({})  -  Jdth(t)e~’Qt  —  e~2Q 
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tion  to  the  target-plane  coherence  length  pi  =  X0L/7ra0.  This  is  the  same  spatial  resolution 
that  was  previously  found  via  Gaussian-# tate  analysis  for  far-held  transmissive  ghost  imaging 
with  a  pseudothermal  source  [4].  Indeed,  the  only  difference  between  Eq  (2.15)  and  the  cor¬ 
responding  result  for  the  transmissive  case  is  the  factor  A2jL2  that  appears  in  the  former. 
In  transmissive  ghost  imaging  all  the  light  that  passes  through  the  target  is  collected  by 
the  bucket  detector,  but  the  quasi- Lambertian  nature  of  the  rough- surfaced  target  combines 
with  the  stand-off  measurement  by  the  bucket  detector  to  introduce  the  solid-angle  subtense 
factor  A2/L2  1  in  Eq  (2.15). 


Contrast  Turning  now  to  the  image  contrast  implied  by  Eq  (2.13).  we  will  employ  the 
contrast  definition  from  [4],  viz. 

c  =  maxK,[C(pi)]  -  min-fc[C(pi)]  ^ 

Co 

with  the  assumption  that  the  target  is  entirely  contained  within  a  region  1Z  centered  at  the 
origin  in  transverse  coordinates  and  having  a  diameter  that  is  much  smaller  than  aL.  For 
simplicity,  we  will  also  assume  that  pL  is  small  enough  to  resolve  all  features  in  the  target's 
intensity-reflection  coefficient,  so  that 

J dp2T(p2)e- «  7 xplTipx),  (2.17) 

and  we  will  take  maxTC[T(pi)]  =  1,  miipj[T(pi)]  =  0.  Thus 

C  «  7T p2L/Ar,  (2.18) 

where  At  =  f  dp>T(p2  )  is  the  effective  area  of  the  target.  Thus  C  ss  1/number  of  011- 
target  resolution  cells.  This  image  contrast  coincides  with  what  was  previously  derived  for 
DC-coupled  transmissive  ghost  imaging  in  far- field  operation  in  [4]. 
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2.3  Signal-to-Noise  Ratio 


The  featureless  background  that  we  encountered  in  the  preceding  section  can  be  eliminated 
by  means  of  AC-coupling  one  or  both  of  the  photodetectors  in  the  Fig.  2-1  setup,  as  has 
sometimes  been  done  in  pseudothermal  ghost-imaging  experiments  [16].  SNR  analysis  for 
transmissive  ghost  imaging  is  simplified  substantially  by  inclusion  of  such  AC-coupled  detec¬ 
tors  [6],  so  we  shall  take  the  same  route  here  by  assuming  that  the  phot.odetectors:  baseband 
frequency  response  include  a  DC  block.  Specifically,  the  detectors’  frequency  response  is  now 
modeled  by 

HB(tt)  =  e~2Q2/n*  -  e“2n2/n^’  (2.19) 


where  the  DC- notch  bandwidth,  0,v.  is  much  smaller  than  both  QB,  the  detectors’  high- 
frequency  cutoff,  and  1/To,  the  source  bandwidth.4  With  this  AC-coupling  we  have  that  the 
average  photocurrent  cross  correlation  is  background  free.  viz..  Eq  (2.13)  becomes 


(G(p  i)> 


q2n2AiA2 

L2 


J  dp2T(P2)e-lpi-p2]2/p c 


(2.20) 


Thus  it  is  appropriate  to  define  the  ghost  image’s  SNR  at  the  image  point  pi  via 


(c{Pl)f  (c(pi)>2 

Varies,)]  (CHPl))  -  (C(Pl)r 


(2.21) 


i.e,  it  is  the  ratio  of  the  squared  strength  of  the  image  component  of  the  photocurrent  cross 
correlation  divided  by  the  variance  of  that  cross  correlation. 

Equation  (2.20)  provides  an  expression  for  the  numerator  in  Eq  (2.21)  and  the  second  term 
in  its  denominator.  However,  to  simplify  our  results,  we  shall  assume  that  pL  is  sufficiently 
small  to  resolve  all  features  in  Tip),  reducing  Eq  (2.20)  to 


(C(Pi)) 


AlAoqY-  (  2 P 


L2 


na. 


KPlTiPi). 


(2.22) 


4Because  we  have  assumed  a  narrowband  pseudothermal  source,  f I^Tq  <C  1  is  a  more  stringent  condition 
than  nN  <C  fffi. 
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This  leaves  us  with  the  formidable  task  of  evaluating 

n'JT'j  j*}  'J'j  jO 

(C2(pi))  =  ~[  dt  f  du  {ii(t)i2(t)ii{u)i2{u)) ,  (2.23) 

1  I  J-T,/ 2  J-T,/ 2 

which  requires  us  to  determine  an  eighth  moment  of  the  fields  and  a  fourth  moment  of  the 
target's  field-reflection  coefficient.  Fortunately,  Gaussian-moment  factoring  can  be  applied 
to  both  of  these  moment  evaluations,  but  the  Fresnel- propagation  kernels  that  canceled  out 
in  finding  the  average  photocurrent  cross  correlation  do  not  do  so  here.  We  can  simplify  the 
analysis  by  using  the  dimensionless  difference  coordinate  v  —  PlA{p>  —  p”)/L,  where  p‘  and 
p"  are  coordinates  at  the  bucket  detector,  and  defining 

4  =  JdPr2(p).  (2.24) 

4'  =  JdHm)\20(SAa/pL).,  (2.25) 

T  =  / civ  e~\A~/20(vAa).  (2.26) 

2ff  J 

where  O(C-  D)  is  the  dimensionless  version  of  the  two-circle  overlap  function  for  circles  of 
diameter  D. 


O(C.D) 


(2.27) 


Q  =  v47  /rag.  and  T(£)  =  J-[T(p)}(^)  is  the  Fourier  transform  of  the  target.  T  was 
normalized  so  that  it  approaches  one  for  very  large  receiving  apertures;  i.e.,  lima  >O0  F  =  1. 
Next,  we  assume  that  A'r/ p\  >  30,  which  is  equivalent  to  saying  that  the  ghost  image 
consists  of  at  least  10  x  10  resolution  cells.  The  full  derivation  under  these  conditions  is 
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shown  in  the  Appendix,  with  the  resulting  SNR  being 


SNR  = 


a't+a'f  ,l7r°2 

%/2 irp? 


T>l) 


Tt  I'^ran 
T0  2Ao 


T\pi) 


Tj_ 

To 


T(p i)  L* 
t;X  a2 


■xpITHpi) 

A\r)X 


4 

3 


a 


2 

0. 


4x^2 


+ 


T(pi)T0nBpj  y^F  J_,2 

H'.v2,l,/pr'  4j 


(2.28) 


where  I  =  PT0p2L/a2L  is  the  source  brightness  in  photons  per  spatiotemporal  mode. 

Equation  (2.28)  was  written  in  its  full  form  to  allow  for  any  size  bucket  detector.  This 
is  useful  when  looking  at  small  bucket  detectors,  such  that  A2  — >  0;  even  for  high-brightness 
illumination  and  a  long  integration  time,  we  are  limited  in  such  situations  to  SNR  <  1 
because  there  is  no  averaging  of  the  target-induced  speckle.  However,  for  a  lensed  ghost 
imaging  system  it  is  reasonable  to  assume  that,  the  area  of  the  collecting  lens  in  front  of  the 
bucket  detector  is  at.  least  as  large  as  the  area  of  the  source  beam,  or  A2j ttOq  >  1.  In  this 
regime,  we  can  simplify  Eq  (2.28)  to 


SNR 


T2(Pi) 


f >  r.T"2 


ppAj _ 

T(pi)  L-  ,  W,.T*(Pi) 
pX  7E  3/li  pI 


TMWbpIV^  Ly 
1  gV2AiV2X2  ^2 

(2.29) 


The  terms  in  the  noise  denominator  of  Eq  (2.29),  which  originate  from  different,  combina¬ 
tions  of  field  variations  and  shot  noises,  have  important,  physical  interpretations.  From  left 
to  right  in  that  denominator  we  have:  the  noise  contributed  by  target-plane  speckle  from 
the  pseudot  henna!  illumination;  the  noise  contributed  by  the  speckle  on  the  bucket  detec¬ 
tor  arising  from  the  target’s  surface  roughness;  the  beat  noise  between  the  pseudothermal 
speckle  on  the  CCD  pixel  and  the  bucket  detector’s  shot  noise;  the  beat  noise  between  the 
CCD  pixel’s  shot  noise  and  the  pseudothermal  speckle  on  the  bucket  detector;  and  the  beat 
noise  between  the  shot  noises  on  the  two  detectors.  From  here  it  is  of  interest  to  look  at  the 
low- brightness  (I<  1)  and  high-brightness  (I>  1)  SNR  asymptotes.  These  are  given  by 


SNRl  = 


16 V2T,  Apft2  A2 

^  T0QBT0pl,[Pl)L2: 


(2.30) 
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and 


SNR 


H  — 


r\P  i) 


'Ll 

To 


+  T*(Pi) 


Tj  r  it  a  ft 

T0  2/1 2 


(2.31) 


respectively. 

The  low- brightness  SNR  is  dominated  by  the  beat  noise  between  the  detectors’  shot 
noises.  It  coincides  with  the  low-brightness  SNR  found  for  transmissive  ghost  imaging  [6], 
except  for  the  following  two  differences:  the  reflective  case  has  the  target’s  average  intensity- 
reflection  coefficient,  T(p\),  appearing  in  lieu  of  transmissive  target’s  |T(pt)|2:  and  the 
reflective  case  includes  the  solid-angle  scaling  factor.  A2/L2.  previously  encountered  in  our 
comparison  of  the  these  ghost  imagers’  spatial  resolutions. 

■  The  reflective  ghost  imager’s  high-brightness  SNR  asymptote  is  controlled  by  the  two 
speckle  terms  from  Eq  (2.29),  i.e..  the  speckle  arising  from  the  pseudothermal  source’s  spa¬ 
tial  incoherence  and  the  speckle  arising  from  the  target’s  surface  roughness.  Neither  speckle 
noise  can  be  said  to  universally  dominate  the  high-brightness  SNR  asymptote,  as  their  rel¬ 
ative  strengths  are  governed  by  both  spatial  and  temporal  factors.  We  need  to  look  at  two 
limiting  cases:  when  the  integration  time  is  short  enough  that  the  source's  spatial  incoherence 
dominates  the  noise,  and  when  the  integration  time  is  long  enough  that  the  target-induced 
speckle  dominates  the  noise.  These  short  integration-time  and  long  integration-time,  high- 
brightness  SNR  asymptotes  are 


rr i  2 

SNR//.  short_T,  =  (Pl),  (2.32) 

l  o  at 

and 

SNR//.  iong-7>  =  Yna2  (2.33) 

Here  we  see  the  short  integration-time,  high-brightness  SNR  for  reflective  ghost  imag¬ 
ing  equals  the  high-brightness  SNR  for  transmissive  ghost,  imaging  with  7^(pi)  appearing 
instead  of  |T(pi)|4.  This  agreement  is  to  be  expected,  as  both  of  these  SNRs  are  limited 
by  the  speckle  created  by  the  pseudothermal  illumination.  However,  as  the  integration  time 
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increases,  the  high-brightness  SNR  for  the  reflective  case  saturates  at  the  value  given  by 
Eq  (2.33).  Here  the  SNR  is  limited  by  the  target-induced  speckle.  Because  we  have  assumed 
a  stationary  target  whose  field-reflection  coefficient,  is  constant,  no  amount  of  post-detector 
integration  will  reduce  its  speckle  noise,  and  SNR  saturation  occurs.  Furthermore,  this  effect 
can  be  severe:  for  A2/'Ka20  =  1  we  find  SNRH,longT,  =  3.266,  and  for  A2/na%  =  2  we  have 
SNR//.iongr,  =  5.54.  So.  for  realistic  stand-off  sensing,  the  SNR  will  be  limited  to  single-digit 
values  if  no  further  measures  are  taken  to  average  out  the  target-induced  speckle. 
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Chapter  3 


Ghost  Imaging  with  Non-Classical 
Light 


3.1  Theoretical  Setup 


Our  analysis  of  reflective  ghost  imaging  in  Chapter  2  focused  on  a  pseudothermal  light  source 
with  classical  phase-insensitive  correlation  between  the  signal  and  reference  fields.  This  setup 
was  chosen  because  we  felt  it  most  closely  resembled  a  system  that  could  be  reliably  con¬ 
structed  for  experimentation.  However,  it  behooves  us  to  consider  alternative  light  sources, 
especially  those  with  quantum  entanglement,  between  the  signal  and  reference  fields.  In  this 
chapter  we  derive  the  resolution,  contrast,  and  SNR  for  ghost  imagers  with  phase-sensitive 
correlation  between  the  signal  and  reference  fields.  We  work  within  a  Gaussian-state  frame¬ 
work  as  done  in  [4],  exploring  correlations  ranging  from  the  classical  limit  up  to  the  quantum 
limit.  But  first,  we  need  to  reframe  our  setup  in  terms  of  quantum  mechanics,  as  shown  in 
Figs.  3-l(a),  3-l(b),  and  3-l(c). 

We  start  by  replacing  the  signal  and  reference  fields  at  the  source  with  their  associated 
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(a)  The  output  from  an  SPDC  is  passed  through  a  polarizing  beam  splitter, 
creating  signal  and  reference  fields  from  the  signal  and  idler  photons  which 
have  a  non-classical  cross  correlation. 


(b)  For  sub-unity  quantum  efficiency  detectors  we  (c)  Since  some  light  is  transmitted  through  the  tar- 
have  to  account  for  the  loss,  and  inject  vacuum  to  get,  the  reflected  field  has  suffered  loss.  We  there- 
preserve  the  field  commutator  relationship.  This  re-  fore  inject  vacuum  as  in  Eq  (3.4),  which  is  similar 
lationship  is  given  in  Eq  (3.G),  and  is  analogous  to  to  the  field  and  vacuum  passing  through  a  spatially 
the  field  operator  and  the  vacuum-state  operator  be-  varying  beam  splitter  of  transmissivity  T(p). 
ing  passed  through  a  beam  splitter  of  transmissivity 

\fn- 

Figure  3-1:  The  setup  in  Fig.  2-1  is  reframed  in  terms  of  quantum  mechanics.  The  source 
is  changed,  the  classical  fields  become  field  operators,  and  we  have  to  inject  vacuum  at  the 
detectors  and  target  to  preserve  commutator  brackets. 
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field  operators  Es(p.t )  and  ER(p,t).  These  fields  obey  the  commutator  relations 


(3.1) 

(3-2) 


[Ei(pi,ti),Em(p2,t2)]  =  0, 


[Et(pi,  t),  Ell{p2-.  t2)\  =  5f,mi5(Pi  -  P2)S{h  -  t2), 


where  t  m  E  {S,R,T}.  The  propagation  of  the  classical  fields  in  Eq  (2.3)  becomes  its 
operator  counterpart 


E'm{p-t)^  /  dp'Ern(p'.t) 


kQ  eiko(L+\p-p'  |2/2L) 

i‘hxL 


for  w  =  R,  S,  T. 


(3-3) 


The  signal  and  reference  held  operators  have  the  same  basic  interpretation  as  their  classical 
counterparts,  i.c.,  having  the  statistics  of  zero-mean  Gaussian  random  processes,  but  the 
target-return  held  description  has  changed.  In  order  to  preserve  the  held  commutator  re¬ 
lations  in  Eqs  (3.1)  and  (3.2),  we  have  to  inject  vacuum  to  compensate  for  the  loss  at  the 
target,  as  shown  in  Fig.  3-  1(c)  and  modeled  as 


ET(p.  t )  -  E's(p,  t)T(p)  +  Evac.s(p.  t } \/ 1  ;/>T.  (3.4) 


where  Evac.s{p,t )  is  a  vacuum-state  held  operator.0 

The  photodetectors  can  now  be  thought  of  as  making  a  measurement  of  a  quantum  held 
operator.  Since  we  are  dealing  with  sub-unity  quantum  efficiency  detectors,  we  have  to  inject 
vacuum-state  held  operators  to  maintain  the  held  commutator  relationships  at  the  detectors, 
as  shown  in  Fig.  3-1  (b) .  Thus,  our  photocurrent  is  a,  measurement  of  the  quantum  operator 

ix{t)  —  q  J dr  J  dpEl(p.  r)Ex(p,  r)h(t  —  t)  for  x  =  1.2,  (3.5) 


BStrictly  speaking,  this  expression  requires  |T(p)|2  <  1  for  all  p,  which  conflicts  with  the  statistics  we 
have  assumed  for  the  field  reflection  coefficient.  However,  the  A2/L  angular  subtense  factor  that  we  will 
encounter  on  the  target  return  will  make  our  statistics  a  reasonable  approximation  for  the  detected  field 
operator.  Also  note  that  the  vacuum-state  field  operator  will  not  contribute  to  the  bucket-detector  output. 
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where 


Ez(p-t)  =  y/i )E'm(p,t)  +  v'l  -1 lEvacm(p.t), 


(3.6) 


and  x  =  1  for  m  =  R.  and  x  =  2  for  m  —  T.  For  this  treatment  our  source  will  be  an 
SPDC,  whose  output  is  Gaussian-state  light  with  a  phase-sensitive  cross  correlation  between 
the  signal  and  reference  (signal  and  idler)  fields  given  by  [5,  6] 


9  P 

(ER(pi,t1)Es{p2,t2))  =  —  e-(led2+|pT)A>§ 

TTOo 

x  |  e-\p2-pi\2/2pjte~{t2-t])V/2T$  ,■  (  ‘^9 

7 T 


PToP20 


h _ _e~\P2—P\  |2/Pog-(i2-0  )2/Tq 


(3.7) 


The  fields  will  be  taken  to  have  no  phase-insensitive  cross  correlation,  no  phase-sensitive 
autocorrelation,  and  a  phase-insensitive  autocorrelation  given  by 


{Eln{pl-tl)Em{p2R2)) 


_ e-(l/>il2+|p2|2)/og-|pi-p2|2/2p§e-(t2-^i)2y^  m  =  R  S  (2  8) 

TTog 


Looking  at  the  Gaussian-state  cross  correlation  function  in  Eq  (3.7)  we  see  two  separate 
correlations.  In  the  low-brightness  limit  the  second  term  becomes  dominant,  and  the  state 
can  be  approximated  as  a  biphoton  state  with  a  correlation  much  stronger  than  the  clas¬ 
sical  limit.  Conversely,  in  the  high-brightness  limit  the  first  term  dominates,  and  the  light 
approaches  a  classical  state,  with  the  correlation  being  the  same  as  for  the  Gaussian- Schell 
model  in  Chapter  2.  Equation  (3.7)  captures  this  full  quantum-to-classical  behavior,  allowing 
us  to  say  that  at  low-brightness  we  have  an  entangled  biphoton  state,  and  at  high-brightness 
we  have  a  classical  state  with  phase-sensitive  cross  correlation  between  signal  and  reference 
fields. 
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3.2  Resolution  and  Contrast 


To  evaluate  the  resolution  and  contrast,  we  follow  t-lie  same  procedure  as  for  classical  light; 
i.e.;  we  look  at  an  ensemble  average  of  the  product  of  the  photocurrents  as 

(C{pi))  =  {tiUiioil)) 

=  q1  A\  J dr}  J d,Toh(t  -  Ti)h(t  -  t2)  J  dp'  {E[{pi,  r1)£’](p'.  r2)Ei{pi,  Ti)E2{p' .  t2)). 

Substituting  in  Eq  (3.6)  and  recognizing  the  independence  of  the  zero-mean  vacuum  fluctu¬ 
ations  from  the  held  fluctuations,  we  have 


{E\(pi,  7 i)Et(p',  r2)£i(pi •  Ti)E2{p' ,  r2))  —  rj~(hR(pi,  r{)Ej-{p  .  T2)ER(p\ ,  T\)ET{p  ,  t2)). 


From  Rq  (3.3)  we  then  find  that. 


{Er{Pi,  h)£(i(Prri)4(p'  ^2)) 


E2 


4ir2E- 


ik.o{\p'-p3 


\P'-P2\2)/2L 


X  ( E R  ( p  1 ,  n )  i?2  ( p2 .  r2 )  E'r  ( p  1 ,  n ) . Eg  ( p 3 ,  To ) )  ( T *  ( p 2 )  T  ( p 3 ) )  ■ 


once  again  leaving  us  needing  to  evaluate  the  fourth  order  held  moment  after  L-meters  ol 
propagation,  as  well  as  the  second  order  target-surface  moment.  The  target  moment  is 
treated  in  the  same  manner  as  it  was  for  classical  illumination,  and  is  thus  evaluated  with 
Eq  (2.5).  Since  our  fields  are  normally  ordered,  and  we  have  Gaussian-state  light,  we  apply 
the  Gaussian-moment  factoring  theorem  to  arrive  at 


(Er(pi,  Ti)ER(p,  T2)E'R(pi.  Ti)E'T{p  ,t2)) 


Ao 

7? 


dp2T(p2 )  (Wr(P\-,e)\2) {\E's(p2,T2)\~)  +  \  (E'R(p1,  Ti) E's(p2.  t2))\ 


which  is  similar  to  the  form  we  found  for  pseudothermal  light,  except  that  we  have  a  phase- 


sensitive  cross  correlation  in  lieu  of  a  phase-insensitive  cross  correlation.  To  finish  the  eval- 
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uation,  we  need  to  propagate  the  phase-sensitive  cross  correlation  and  phase  insensitive 
autocorrelation  L  meters  as  follows: 


(E'R{pi,tl)E's(p2.:t2)) 


( e£\pi\2e£\p*\2e-\K+p?\2/24 

\™l ) 


X 


e-(\m?+M2)/«2oe^-t,)2/2Tg  +  1 


(|Pi!2  +  iP2j2)/2a5e 


{t2-U)2/t$ 


(3.9) 


and 


{E'l{piJ.1)E'm(p2,t2)) 


2  P 


1KQ 

-e-2L 


(|P2i2-|pi|2)p-(|pi|2  +  [p2!2)/a'i-|pi-p2|2/2p|  -(i2-O)2/2r02 


na 


(3.10) 


for  m  =  R.  S.  We  can  now  evaluate  the  ensemble  average  of  the  correlation  function,  once 
again  assuming  that  the  on-target  average  illumination  pattern  is  wide  enough  that  it  is 
essentially  uniform  across  the  target.  With  that,  we  find  that  the  mean  of  the  correlation 
becomes 


(CM) 


r/-y.b  .1(1  (  2 P 


L- 


7 ra~, 


dpi  Tip  , 


1  +  e"#1  +pT/pI  ( 1 _ l  1 

\2yfal  , 


(3.11) 


where  once  again  J  =  PT0p2L/a2L.  As  was  done  in  Chapter  2.  this  can  be  broken  into  the 
background  term 


Ch  = 


qWAiA2  f  2 P 


L2 


ncr. 


dpi  Tipi). 


(3.12) 


and  the  image-bearing  term 


CApx)  =  —  (K)  f  dp,  T(p,)e-'^fd‘ln  (  1 


l2  V™!, 


+1 


(3.13) 


Resolution  Comparing  this  to  our  results  from  pseudothermal  illumination,  we  see  that 
the  image  is  blurred  by  the  same  PSF,  except  that  the  coordinates  are  inverted  (pi  — »  —pi). 
Thus,  the  image  has  the  same  resolution  of  p^  =  X^L/nao. 
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Contrast  Comparing  Eq  (3.12)  to  Eq  (2.14),  we  see  that  phase-sensitive  and  phase- 
insensitive  imaging  have  the  same  background  term.  Comparing  Eq  (3.13)  to  (2.15),  we 
see  that  the  image-bearing  terms  differ;  this  is  due  to  our  inclusion  ot  non-classical  correla¬ 
tions  effects.  If  we  look  at  the  high-brightncss  case,  the  classical  correlation  dominates  the 
non-classical  correlation,  and  the  image-bearing  term  for  classical  phase- sensitive  imaging  is 
the  same  as  the  image-bearing  term  for  phase-insensitive  imaging  (except  for  the  coordinate 
inversion). 

To  evaluate  at  the  contrast  we  will  use  Eq  (2.16),  the  definition  we  used  for  the  phase- 
insensitive  imaging.  Using  the  simplifying  assumption  that  the  PSF  is  narrow  enough  to 
resolve  all  target-features,  and  that  maxK[T(pi)]  =  1  and  niin-R.[T(pi)]  =  0,  we  can  say 

J df>-2Ti(>2)<  ys  rrpyTU-p:).  (3.14) 


and  simplify  our  contrast  definition  to 


C 


*Pl 


(3.15) 


where  once  again  AT  =  Jdp2T(p2)-  In  the  high-brightness  limit,  Eq  (3.7)  approaches  a 
classical  phase-sensitive  cross  correlation  and  the  contrast  in  Eq  (3.15)  approaches  Eq  (2.18). 
our  result  from  the  pseudothermal  ghost  imager.  However,  in  the  low-brightness  limit  the 
the  quantum  correlation  dominates,  and  the  contrast  in  Eq  (3.15)  grows  without  bound. 
This  is  the  same  behavior  seen  in  transmissive  biphoton  imaging  [1,  4].  Thus,  maximally 
correlated  phase-sensitive  imaging  has  the  same  resolution  as  phase-insensitive  imaging,  but 
improved  contrast  in  the  low-brightness  limit. 


3.3  Signal-to-Noise  Ratio 

The  derivation  of  the  SNR  for  phase-sensitive  light  will  follow  what  was  done  in  Section  2.3 
for  the  phase-insensitive  case.  The  SNR  will  be  defined  as  in  Eq  (2.21)  as  the  ratio  of  the 


33 


squared  mean  to  the  variance  for  the  image-bearing  portion  of  the  photocurrent  correlation, 
which  will  now  be  given  by  Eq  (3.13).  We  will  be  implementing  a  DC-block  in  the  detector 
response  as  modeled  by  Eq  (2.19),  and  will  assume  that  the  target  features  are  completely 
resolved  by  the  imager,  allowing  us  to  use  the  simplification  in  Eq  (3.14).  With  that,  the 
mean  of  the  image- bearing  term  becomes 


Ci  ( Pi  | 


q2rj2AlA2 

L2 


nplT(~Pi) 


leaving  us  needing  to  evaluate  the  image-bearing  term's  second  moment 


l  rTr/2  rTf/  2 

(C2(pi))  =  7^  /  dt  du(z'1(t)i2(t)ii(u)?2(n)). 

LI  J-T,/ 2  7-77/2 


(3.16) 


(3.17) 


In  evaluating  Eq  (3.17),  wc  have  our  first  major  departure  from  the  derivation  in  Section 
2.3.  In  Chapter  2  we  used  a  semiclassieal  treatment  of  the  light,  where  fluctuations  in 
the  measurement  were  treated  as  coming  from  the  conversion  of  the  continuous  field  to 
discrete  charges  in  the  detector,  which  we  called  the  detector  shot  noise.  However,  we 
are  now  dealing  with  quantum-mechanical  entangled  fields,  and  must  use  a  full  quantum 
treatment.  In  quantum  optics  theory,  measurement  fluctuations  arise  from  fluctuations  of 
the  quantum  field,  coupled  with  the  type  of  measurement  being  performed.  Mathematically, 
these  noise  terms  appear  from  applying  the  commutator  operations  in  Eqs  (3.1)  and  (3.2) 
when  normally  ordering  our  higher  order  field- moments.  While  the  noise  qualitatively  comes 
from  different  sources,  quantitatively  this  has  the  same  effect  as  the  semiclassieal  treatment: 
after  normally-ordering  our  fields  we  have  an  eighth-order  moment,  two  sixth-order  moments, 
and  a  fourth-order  moment.  This  comes  as  no  surprise,  as  it  was  shown  in  [4]  that  ghost- 
image  formation  using  a  downconversion  source  is  inherently  a  classical  phenomenon,  with 
the  only  non-classical  features  coming  from  the  stronger-than-classical  correlation  of  signal 
and  idler  photons.  As  before,  since  these  higher-order  moments  are  all  normally  ordered,  and 
the  fields  are  zero-mean,  we  can  apply  the  Gaussian-moment  factoring  theorem  to  express 
them  as  the  sum  of  products  of  second  order  moments. 
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We  evaluate  the  moments  at  a  distance  of  L-meters  from  the  source,  using  Eq  (3.9)  for 
the  phase-sensitive  cross  correlation  and  Eq  (3.10)  for  the  phase-insensitive  autocorrelations. 
All  other  second-order  moments  evaluate  to  zero.  We  will  again  move  to  the  normalized 
and  difference  coordinate  v  =  pzko(p'  —  p")/L ,  where  p'  and  p"  are  coordinates  at  the 
bucket  detector.  The  definitions  of  T,  a,  A'F,  and  A'T  are  the  same  as  in  Chapter  2.  Finally, 
assuming  that  there  are  at  least  10  x  10  resolution  cells  in  the  image  and  inverting  our 
coordinate  system  as  px  — ¥  —p\,  we  have 


SNR  = 


4'  |  "oa'f 

4ttAo 


v^p2lT(piP 


T,  meg  p 

T02A2T(pt)2 


-I2 


Kpj 

MTrjT(pi)  1  A\2r) 


L2 


+ 


4ttA2 


•/2L2T0nDpl 


I 


16  A2AiA2Tir)2T\pi)) 

(3.18) 


where  I  =  (777^5  +  *•)  ■  Please  see  the  Appendix  for  the  full  derivation.  This  result  is  quite 
unwieldy,  so  we  would  like  to  simplify  it  by  assuming  that  the  receiving  aperture  is  at  least 
as  large  as  the  source  aperture,  giving  us 


SNR 


A'r  ,  TiPl)L-  W;r-(p,)  T(pi)y^7o.QBpi  ry 

[pi  T0  2A2  ip  A2  SA^T)  a2 

(3.19) 


This  result  corresponds  to  the  transmissive  result  calculated  in  [6].  in  the  same  manner  as 
its  classical  phase-sensitive  counterpart  in  Chapter  2  did;  that  is,  there  is  now  a  solid-angle 
subtense  factor  on  terms  associated  with  the  bucket  detector,  and  a  target-speckle  term  in 
the  variance  that  does  not  diminish  with  integration  time. 


In  the  classical  limit  I  — >•  1,  Eq  (3.19)  becomes  Eq  (2.29),  and  we  can  see  that;  classically 
correlated  phase-sensitive  imaging  has  the  same  SNR  characteristics  as  phase-insensitive 
imaging  with  pseudothermal  light.  This  means  that  the  high-brightness  limit  is  also  the 
same,  and  is  given  by  Eq  (2.31).  Now,  in  the  low-brightness  limit  the  quantum  correlation 
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significant.  I  -4  and  the  low-brightness  SNR  becomes 

Ti  8 AirfX  M 

t0mbt0Pi  MPlV 


(3.20) 


This  equals  the  classical  low-brightness  limit  divided  by  a  factor  of  2\/27rZ,  which  means 
that  for  I<1,  the  low-brightness  SNR  is  orders  of  magnitude  larger  for  quantum-correlated 
light  than  for  classically  correlated  light.  This  is  the  second  signature  of  biphoton  state  ghost 
imaging:  improved  low-brightness  SNR. 
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Chapter  4 


Computational  Reflective  Ghost 
Imaging 


4.1  Theoretical  Setup 


Ghost  imaging  requires  knowledge  of  the  time-varying  speckle  pattern  illuminating  the  tar¬ 
get.  Because  the  ground-glass  diffuser  in  Fig.  2-1  randomly  modulates  the  source  field,  we 
measure  the  reference  field’s  speckle  pattern  with  the  CCD  array,  and  exploit  its  correlation 
with  the  speckle  pattern  impinging  on  the  target  to  form  the  ghost  image.  Suppose,  however, 
that  a  known  source  is  subjected  to  a  deterministic  spatiotemporal  modulation,  through  use 
of  a  spatial  light  modulator  (SLM),  in  a  manner  that  projects  a  time- varying  but  determin¬ 
istic  speckle  pattern  on  the  target.  In  this  case  the  speckle  pattern  at  the  target  can  be 
computed  from  diffraction  theory,  and  we  do  not  need  the  reference  arm  to  form  a  ghost 
image  [9],  As  shown  in  Fig.  4-1,  we  form  a  computational  ghost  image  by  cross  correlating 
the  computed  reference-arm  photocurrent  , 

c(£)  =  J 


dr  qr]Pi(T)h(t  —  r), 


(4.1) 


Bucket 


Figure  4-1:  Setup  for  computational  reflective  ghost  imaging. 


with  the  measured  photocurrent,  hit),  from  the  bucket  detector.  Here. 


Pi(t)  =  AxlEiipiJ)]2, 


(4.2) 


gives  the  computed  photon  flux  for  a  CCD  pixel  located  at  pi  in  terms  of  the  computed 
speckle  pattern  |  /-  i  ipi  .  C|'. 


The  SLM  is  traditionally  a  square  of  size  D  x  D,  composed  of  a  grid  of  square  pixels 
of  size  d  x  d.  Each  pixel  centered  at  transverse  location  pt  causes  a  phase  shift  <pt{t)  on 
the  light  impinging  it,  imparting  a  spatially- varying  phase  shift  on  the  resulting  field.  This 
means  that,  assuming  the  light  hitting  the  SLM  can  be  approximated  as  a  plane-wave,  the 
field  leaving  the  SLM  is 


E(p.t)  =  ]T>ct  reel  i^Aj 


(4-3) 


where  rect(-)  is  the  unit-lengtli  rectangle  function,  p  =  (x,  y),  and  pt  =  [xp  yi)  is  the  center 
of  pixel  i.  Propagating  this  into  the  far  field  results  in 


d2 


x'd 


y'd 


E'(p',t)  =  V  -——sine  (  )  sine  (  )  e-^(x'xe+y'ye)/Le-zk0(\x'\2+\y'\ 2)pL  j<pt{t) .  (4.4J 

'  A oL  V  A f)L  1 


An  L 


where  sinc(-)  is  the  sine  function,  and  p'  —  (x'.y'). 
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4.2  Comparison  to  Pseudothermal  Ghost  Imager 


For  comparison  with  the  pseudothermal  ghost  imager,  we  would  like  to  derive  the  spatial 
resolution,  image  contrast,  and  SNR  of  this  computational  counterpart.  Looking  at  Eq 
(4.4),  it  is  not  immediately  clear  how  to  make  the  connection  to  imaging  with  Gaussian- 
state  light.  However,  as  the  field  in  Eq  (4.4)  is  the  sum  of  a  sufficiently  large  number 
of  weighted  independent,  identically  distributed  (IID)  random  variables,  from  the  Central 
Limit  Theorem  it  will  have  Gaussian  random  process  statistics  [9].  This  can  be  seen  in 
Fig.  4-2,  which  shows  a  pseudorandom  phase  modulation  pattern  in  (a),  and  the  resulting 
far- field  intensity  pattern  in  (b).  Thus,  we  can  treat  the  fax-field  pattern  E's(p.  t)  hitting  the 
target  and  the  reference  field  Ei(pi,t)  =  E'R(p,  t)  —  now  calculated  from  the  phase  shifts 
by  Eq  (4.4)  —  as  zero-mean,  complex- valued  Gaussian  random  processes.  Furthermore,  to 
simplify  our  analysis,  we  shall  assume  that  we  can  use  the  Gaussian-Schell  model  for  the  field 
correlations  as  given  in  Eq  (2.10).  Under  these  conditions,  the  far-field  coherence  length  can 
be  approximated  as  pL  «  ‘2L/k0D ,  and  the  far-field  intensity  radius  by  aL  «  2 L/kd.  With 
these  assumptions,  all  the  derivations  from  Chapter  2  carry  over  to  computational  ghost 
imaging  by  simple  omission  of  the  CCD  array's  shot  noise. 

Resolution  Under  the  Gaussian-Schell  assumption  the  spatial  resolution  of  the  computa¬ 
tional  ghost  imager  is  identical  to  that  of  the  pseudothermal  ghost-  imager  at  pL  —  2 L/k0d; 
therefore,  in  practice,  the  resolutions  should  be  similar. 


Contrast  The  computational  imager  has  the  same  image  contrast  as  pseudothermal  ghost- 
imaging  for  DC-coupled  operation,  found  in  Eq  (2.18).  Moreover,  AC-coupling  of  the  bucket, 
detector’s  photocurrent  or  the  CCD  array’s  computed  photocurrent  wall  eliminate  the  fea¬ 
tureless  background  term  in  the  photocurrent  cross  correlation,  giving  the  computational 
ghost  imager  the  same  high-contrast,  behavior  seen  earlier  for  AC-coupling  in  pseudothermal 
ghost  imaging. 
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Signal-to-noise  ratio  There  is  an  interesting  difference  between  computational  and  pseu- 
dothermal  ghost  imaging  that  appears  when  we  compare  their  SNR  formulas.  Since  the 
reference  arm  is  computed,  there  are  no  shot  noise  fluctuations  on  the  current  R(£)  asso¬ 
ciated  with  our  high-spatial  resolution  reference  arm,  and  the  computational  ghost  image's 
SNR  is  given  by 


c  j\m 
GIN  ltComp 
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V2ttp2i 


T\p  i) 
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(4.5) 


Comparing  this  formula  with  the  pseudothcrmal  result  from  Eq  (2.29)  we  see  that  the  noise 
denominator  for  the  computational  ghost  image's  SNR  contains,  from  left  to  right,  terms 
that  represent:  the  noise  from  the  speckle  pattern  cast  on  the  target;  the  noise  from  the 
speckle  on  the  bucket  detector  arising  from  the  target’s  surface  roughness;  and  the  beat  noise 
between  the  computed  field  speckle  on  each  pixel  and  the  bucket  detector's  shot  noise.  It 
follows  that  the  high-brightness  SNR  asymptote  for  computational  ghost  imaging  is  identical 
to  that  for  pseudothermal  ghost  imaging,  as  it  is  limited  by  the  source  and  target-induced 
speckle.  However,  the  computational  ghost  image’s  SNR  enjoys  a  considerable  advantage  at 
low  source  brightness,  viz.,  its  low-brightness  asymptote  of 


SNRcomp.  L  =  TipA-lp-r/ljA.  (4.6) 

i  o  L 

which  is  significantly  higher  than  that  for  the  pseudothermal  ghost  imager  because  it  scales 
linearly,  rather  than  quadratically,  with  source  brightness. 

Interestingly,  the  computational  low-brightness  limit  also  compares  favorably  to  the  non- 
classical  limit  in  Eq  (3.20).  Both  scale  linearly  with  the  brightness,  but  given  our  receiver 
bandwidth  and  pixel  size  assumptions  of  OBT0  1  and  A±  -C  p\,  the  SNR  for  computational 
imaging  in  the  low-brightness  limit  is  orders  of  magnitude  larger  than  the  corresponding  low- 
brightness  (biphoton)  limit  for  non-classical  ghost  imaging. 
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(a)  The  pseudorandom  phase  modulations  applied 
to  the  pixels  of  the  SLM.  The  phases  at  each  pixel 
are  IID  uniform  random  variables  on  0  to  2tt,  with 
a  new  realization  at  each  time  epoch. 


Figure  4-2:  One  realization  of  a  pseudorandon 
calculated  far-field  intensity  pattern. 


(b)  The  far-field  intensity  pattern  cast  by  the  SLM. 
The  field  distribut  ion  closely  mimics  that  of  a  Gaus¬ 
sian  random  process,  with  the  coherence  length  be¬ 
ing  inversely  proportional  to  the  extent  of  the  SLM, 
and  the  intensity  radius  inversely  proportional  to 
the  width  of  each  pixel 

1  phase  pattern  cast  on  an  SLM  and  resulting 
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Chapter  5 


Comparison  to  a  Laser  Radar  System 

5.1  Theoretical  Setup 

The  importance  of  ghost  imaging  for  stand-off  sensing  rests  on  it  offering  some  advantage 
over  a  comparable  laser  radar  system  for  the  same  application.  We  will  use  the  results  from 
Chapters  2,  3,  and  4  to  provide  a  preliminary  comparison  between  ghost  imaging  and  laser 
radar  for  stand-off  sensing.  The  laser  radar  system  we  shall  consider  is  shown  in  Fig.  5-1.  It 
is  a  direct-detection  system  in  which  a  cw  laser  beam  is  used  t.o  produce  a  spatially  coherent 
beam  at  range  L- m  whose  deterministic  intensity  pattern  matches  in  both  photon  flux  and 
intensity  radius  the  average  intensity  pattern  of  the  Gaussian- Schell  model  we  employed 
for  the  ghost  imager.  A  fraction  of  the  laser  light  reflected  by  the  target  is  focused  onto  a 
CCD  array  by  a  lens  that  is  co-located  with  the  laser  transmitter.  A  target  image  is  then 
formed  by  Tj-s  time  averaging  of  the  output  currents  from  each  CCD  pixel.  The  entrance 
pupil  for  the  laser  radar’s  receiving  lens  will  be  taken  to  coincide  with  the  bucket  detector's 
Ao  active  region  in  the  ghost  imager,  and  we  will  assume  shot-noise  limited  CCD  operation 
as  was  the  case  for  the  Fig.  2-1  setup.6 


6 We  recognize  that  most  laser  radar  systems  employ  pulsed  sources.  We  have  chosen  the  cw  case  to  put 
the  laser  radar  on  the  most  equal  footing  with  the  ghost  imager  for  a  baseline  comparison  between  their 
spatial  resolutions,  image  contrasts,  and  signal-to-noise  ratios. 
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The  laser  radar  image  for  the  CCD  pixel  at  location  p\  is 

I  rTi/% 

I{pi)  =  7^  dti3(t ),  (5.1) 

I I  J-Ti/2 

where  i3(t)  is  the  photocurrent  from  that  pixel.  This  photocurrent  will  have  the  same 
structure  as  seen  in  Eq  (2.2)  for  ghost  imaging,  i.e.,  it  will  consist  of  a  term  driven  by  the 
photon  flux  P3(t)  illuminating  the  pixel  in  question  plus  the  shot  noise  from  that  pixel.  We 
shall  assume  the  DC-coupled  photodetector  model  from  Eq  (2.11),  and  we  will  assume  the 
pixels  are  small  enough  that  they  do  not  limit  the  laser  radar’s  spatial  resolution.  Once 
again  we  shall  assume  1:1  imaging,  although  the  actual  system  will  cast  a  minified  image  on 
the  CCD,  and  wc  shall  invert  the  image  plane  coordinates  so  as  to  obtain  an  erect  image 
of  the  target.  With  these  assumptions  standard  Fourier-optics  thin  lens  theory  leads  to  the 
following  expression  for  E3(pi,t),  the  field  illuminating  the  CCD  pixel  at  p\ : 


E3{pi,t) 


JdpT(p)e^2/2L 


\p\'2/a2L  ^2  Jl(7rZ)o|p  -  Pi\/X0L) 

A5T2  ttD2\p  —  pi\/2X0L 


(5.2) 


where  Do  is  the  diameter  of  A2,  Ji  is  the  first-order  Bessel  function  of  the  first  kind,  and  we 
have  suppressed  absolute  and  quadratic  phase  factors  that  do  not  contribute  to  \E3(pi,t)\2. 
The  photon  flux  for  the  pixel  at  p1  is  thus  P3(t)  =  Ai\E3(p3,t) |2.  and  the  photocurrent  i3(t) 
is 

c(t)  =  qiiD’M)  +  A  i3(t),  (5.3) 

where  the  second  term  is  the  shot  noise.  Since  the  source  is  not  fluctuating  the  detector 
response  will  not  have  a  noticeable  effect  on  the  photocurrent,  and  has  been  omitted.  In 
keeping  with  what  we  did  for  ghost  imaging,  we  shall  assume  that  the  target  is  uniformly 
illuminated  by  the  laser  radar,  so  that  we  can  use  e_lp|2/«i  ~  1  m  gq  (5.2). 
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space 

propagation 

Figure  5-1:  Setup  for  direct-detection  laser  radar. 


5.2  Spatial  Resolution  and  Image  Contrast 


To  derive  the  laser  radar’s  spatial  resolution  and  image  contrast,  we  once  again  look  at  the 
ensemble-averaged  image.  Averaging  Eq  (5.1)  over  the  target’s  surface  roughness  and  the 
detector’s  shot  noise  we  find 


(/(Pi)> 


qyAlAl  2 P  f  (J^D2\p-Pl\/\0L)\2 

A lLA  7T a2LJ  P  P\  irD‘>\p  Pi|/2A0T  / 


(5.4) 


which  shows  that  the  target  image  is  proportional  to  the  target’s  intensity-reflection  coeffi¬ 
cient  convolved  with  the  familiar  Airy  disk  PSF  for  incoherent  imaging. 


Resolution  Our  laser  radar  has  a  spatial  resolution  given  by  \.22\)L/ D2.  In  our  ghost 
imaging  setup,  with  D2  =  2a0,  the  spatial  resolution  is  given  by  pi  =  A0L/7ra0  -  2\)L//nD2 , 
which  is  comparable  to  that  of  the  laser  radar. 


Contrast  Equation  (5.4)  also  shows  that  our  direct-detection  laser  radar’s  image  is  not 
embedded  in  a  featureless  background,  making  its  image  contrast  superior  to  that  of  DC- 
coupled  pseudothermal  ghost  imaging  but  equivalent  to  that  of  the  AC-coupled  version. 
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5.3  Signal-to-Noise  Ratio 


We  define  the  laser  radar’s  SNR 


SNR  = 


by 

(/(Pi))2 

Varf  I  (p^] 


(I(Pi))2 

(HPi))  ~  {/(Pi))2 


(5.5) 


to  enable  a  direct  comparison  with  the  ghost  imaging  SNR  from  Eq  (2.21).  We  will  assume 
that  the  Airy  disk  PSF  resolves  all  significant  features  in  the  target's  intensity-reflection 
coefficient,  yielding 


,rf„w  QVAiAo2P 

«?,)>  = -y-^T(Pl). 


7ra| 


(5.6) 


The  variance  calculation  we  need  is  much  simpler  than  what  we  performed  for  ghost  imaging. 
Making  use  of  the  iterated-expectation  formula, 


Var[/(pi)]  -  £’{p3(i):_2;/2<f<7>/2}[Var( /(px))  I  {-P3(f)  ;  ~Tj/ 2  <  t  <  Tj/ 2})] 

+  Var{p3(t)._T//2<t<2;/2}[£'(^(Pi)  |  {P'3 (t)  :  —Tjj 2  <  t  <  Tj/2})].  (5.7) 


we  can  easily  evaluate  the  noise  denominator  in  Eq  (5.5).  The  first  term  on  the  right  in 
Eq  (5.7)  is  due  to  the  target-induced  speckle,  and  is  given  by  {/(pi))2,  and  the  second  term 
on  the  right  in  that  equation  is  due  to  the  shot  noise,  and  is  given  by  q(J(p1))/Ti.  This 
leaves  us  with 


SNR 


T(pi) 


Pi  Pi)  + 


L2  Teal 
r)AxA2  2 FT, 


(5.8) 


From  this  SNR  expression  we  immediately  sec  that  when  both  the  source  brightness 
and  integration  time  are  sufficiently  high,  the  laser  radar’s  SNR  saturates  at  a  maximum 
value  of  unity,  limited  by  the  target-induced  speckle.  With  the  small  CCD  pixels  we  have 
assumed,  ghost  imaging  still  experiences  a  spatial  averaging  of  the  target  speckle  on  the 
bucket  detector,  whereas  no  such  effect  is  available  for  the  laser  radar  system.  Thus  the 
laser  radar’s  performance  is  inferior  to  that  of  the  ghost  imager  when  both  systems  have 
target-speckle  limited  SNRs.  Outside  of  this  limiting  scenario,  the  relationship  between  the 
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two  systems’  SNRs  is  more  complicated,  as  we  will  now  show. 

To  compare  the  stand-off  sensing  SNRs  for  ghost  imaging  and  direct-detection  laser  radar 
in  more  detail,  consider  the  plots  shown  in  Figs.  5-2(a)  and  5-2(b).  Both  figures  assume 
A0  =  1.5  /rm  operation  for  ghost  imaging  and  laser  radar,  with  their  transmitters  having  the 
same  photon  flux.  For  ghost  imaging  we  assume  the  source  parameters  are  a0  =  1  cm  and 
p0  =  O.lS/rr  mm.  The  target  is  assumed  to  be  at  L  =  1  km  range,  with  effective  area  A'r  = 
100  m2.  Thus  aL  =  X0L/np0  =  10  in  implies  that  the  target  illumination  is  nearly  uniform 
on  average,  with  speckle-limited  spatial  resolution  given  by  pL  —  \0L/ti%  =  OTS/ttiii. 
The  CCD  array’s  pixel  area  will  be  taken  to  satisfy  Ax  =  0.1  p|,  and  the  bucket  detector’s 
area  will  be  set  to  A2  =  ti%.  Both  detectors  will  have  ij  =  0.9  quantum  efficiency  with 
bandwidths  obeying  QBT0  —  100.  The  laser  radar’s  transmitter  will  produce  a  spatially 
coherent  Gaussian  intensity  pattern  on  the  target  with  the  same  «/.  value,  and  its  CCD 
array  will  be  identical  to  that  of  the  ghost  imager.  We  are  interested  in  the  SNR  behavior 
of  our  systems  as  a  function  of  source  brightness  and  detector  integration  time:  specifically, 
we  look  at  the  detector  integration-time  dependence  of  the  SNRs  for  a  strong  source,  and 
the  SNR  source-strength  dependence  for  a  long  integration  time. 

Figure  5-2(a)  plots  the  SNRs  for  pseudothermal  ghost  imaging,  non-classical  SPDC  ghost 
imaging,  computational  ghost  imaging,  and  direct-detection  laser  radar  versus  the  normalized 
integration  time.  T//T0,  when  T{p\)  =  1  and  the  transmitter’s  source  brightness  is  J  = 
10°  photons /mode.  All  three  gliost-imaging  systems  show  the  same  high- brightness  behavior, 
and  we  see  a  slight  SNR  advantage  for  laser  radar  operation  when  Tj/Tq  <  104,  with  the 
ghost  imagers  offering  higher  SNRs  when  all  three  systems  approach  their  target-speckle 
limits.  Figure  5-2 (b)  plots  the  four  systems’  SNRs  versus  the  source  brightness,  assuming 
T{p\)  —  1  and  7j/T0  =  107.  Here  we  see  that  computational  ghost  imaging  provides 
the  best  performance,  while  laser  radar  operation  is  the  worst  performer  except  for  X  ~ 
1CT2.  Pseudothermal  ghost  imaging  is  outperformed  by  computational  operation  until  both 
systems’  SNRs  reach  their  common  target-speckle  limit.  At  very  low  brightnesses,  the  SPDC 
imager  trails  only  the  computational  imager,  but  its  performance  quickly  converges  with  that 
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(a)  Signal-to-noise  ratios  versus  normalized  integra-  (b)  Signal-tonoisc  ratios  versus  source  brightness, 
tion  time,  Tj/T0,  for  high  source  brightness  1=1=  PT0po/a,%.  for  long  integration-time  operation 
PTQf$/a%  =  109.  Tj/T0  =  107, 


Figure  5-2:  Signal-to-noise  ratio  comparison  between  pseudothermal  ghost  imaging,  compu¬ 
tational  ghost  imaging,  and  laser  radar  operation.  The  parameter  values  assumed  are  given 
in  the  text. 


of  the  pseudothermal  imager. 

Thus,  for  very  low-brightness  illumination,  the  SNRs  arc  ordered  from  best  to  worst  as 
follows:  computational  ghost  imager,  SPDG  ghost  imager,  laser  radar,  and  pseudothermal 
ghost  imager.  However,  it  should  be  noted  the  the  SNRs  in  this  range  are  of  the  order 
10-8,  which  provides  little  usable  information.  Once  the  brightness  has  increased  to  the 
point  where  we  can  retrieve  some  image  information,  the  SNRs  are  clearly  arranged  from 
best  to  worst:  computational  ghost  imager,  SPDC  and  pseudothermal  ghost  imagers  (which 
have  converged),  and  the  laser  radar.  It  thus  is  reasonable  to  say  that,  for  the  systems  we 
have  compared,  that  computational  ghost  imaging  has  the  best  SNR  behavior,  followed  by 
SPDC  and  pseudothermal  ghost  imaging,  with  the  laser  radar  at  the  back.  However,  for 
short  integration  times  and  a  very  high-brightness  source,  there  is  an  advantage  to  the  laser- 
radar  system.  This  is  similar  to  the  behavior  seen  in  pulsed  laser-radar  systems,  indicating 
that  for  a  definitive  comparison  between  ghost  imaging  and  laser  radars,  we  should  include 
comparisons  to  a  pulsed  laser  radar  systems. 
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Chapter  6 


Evaluation  of  Turbulence  Effects 

6.1  Turbulence  Review 

Our  analysis  in  the  preceding  chapters  had  assumed  all  light  fields  were  propagating  through 
free  space;  that  is,  wc  have  not  taken  into  account  any  effects  the  propagation  medium 
might  have  on  real-world  experiments.  In  many  stand-off  sensing  applications  the  fields  will 
be  propagating  through  the  atmosphere.  The  earth’s  atmosphere  is  comprised  of  a  variety 
of  gases  and  entrained  particulates,  causing  wavelength-dependent  absorption  and  Rayleigh 
scattering  [17].  These  effects  will  cause  attenuation,  but  that  can  be  minimized  with  an 
appropriate  choice  of  the  operating  wavelength.  The  most  deleterious  atmospheric  effects 
arise  from  propagating  in  bad  weather,  i.e.,  though  fog  or  clouds  [17].  Here  the  scattering  is 
so  severe  that  neither  laser  radar  nor  ghost  imaging  can  be  expected  to  provide  useful  standoff 
imaging.  However,  even  in  clear  weather,  the  atmosphere  is  in  constant  flux.  This  random 
mixing  of  the  air  parcels  with  ~  1A'  temperature  fluctuations  create  random  spatiotemporal 
variations  in  the  refractive  index  known  as  atmospheric  turbulence  [17].  For  a  real  world 
analysis,  the  turbulence-induced  effects  on  the  resolution  of  a  remote-sensing  system  must  be 
explored.  Initial  work  has  been  done  to  explore  turbulence  for  transmissive  ghost  imaging 
with  classical  pseudothermal  light  in  [11];  our  results  will  be  for  reflective  imaging  with 
classical  and  non-classical  sources,  computational  imaging,  and  also  provide  a  comparison 
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to  laser  radar  systems. 

To  start,  we  shall  consider  operation  in  good  weather,  at  a  wavelength  for  which  absorp¬ 
tion  and  scattering  effects  are  minimal.  In  such  scenarios,  atmospheric  turbulence  yields 
refractive  index  changes  on  the  order  of  lCT0  [17].  While  these  changes  might  seem  small, 
their  effect  on  light  propagating  over  a  long  distance  is  profound.  For  instances,  two  in-phase 
fields  at  A  =  1.5  /Jin  propagating  through  media  with  a  refractive  index  difference  of  just 
10-6  will  be  out  7 r  rad  of  phase  in  just  .75  m.  For  propagations  over  several  kilometers,  with 
refractive  index  variations  in  both  time  and  space,  the  accumulated  phase  distortions  will  be 
significant.  Moreover,  because  these  phase:  changes  are  spatially  varying,  constructive  and 
destructive  interference  occurs,  i.e.,  initial  phase  fluctuations  lead  to  intensity  fluctuations 
known  as  scintillation  [18]. 

For  our  analysis  we  will  bo  using  the  Kolmogorov  model  for  turbulence.  Kolmogorov 
said  that  for  turbulence  that  caused  the  spatial  variations  of  the  refractive  index  to  occur 
on  distance  scales  between  a  maximum  L0  cs  10  —  100  m  and  a  minimum  —  10_3m, 
the  temperature  fluctuation  spatial-structure  function  followed  a  two-thirds  power  law  [19]. 
That  is,  the  second  moment  of  the  temperature  difference  between  two  points  is  proportional 
to  two-thirds  the  distance  between  the  points  as 

Dtt(p)=((T(Pi)  +  p)-T(p0))2) 

•  -  f  'rl/'P-  (6.1) 

where  DTT  is  the  temperature  structure  function  and  (7f.  is  the  temperature  structure  con¬ 
stant. 

The  temporal  fluctuations  come  from  two  sources:  the  changing  of  the  shape  of  the 
refractive  index  structure  as  the  current  eddies  mix  the  air,  and  the  drifting  of  the  structures 
with  the  average  wind  velocity  [20].  In  normal  conditions,  the  time  evolution  is  dominated 
by  the  latter  mechanism,  and  the  refractive  index  structure  thus  has  a  typical  coherence  time 
of  tc  =  10-3  —  10“2  second.  This  means  that  on  shorter  time  scales  the  turbulence  can  be 
thought  of  as  frozen. 
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There  has  been  a  lot  of  work  done  in  analyzing  laser  beam  propagation  through  turbu¬ 
lence,  and  for  our  work  we  take  advantage  of  the  Extended  Huygens- Fresnel  Principle  [17]. 
This  is  written  as  the  propagation  of  a  spherical  wave,  multiplied  by  some  path-dependent 
complex  exponential  which  encompasses  the  turbulence-induced  amplitude  and  phase  vari¬ 
ations  for  a  particular  path,  defined  by  its  starting  and  ending  coordinates.  For  Xm{p-  p') 
and  (prn  ( p •  p')  being  the  log- amplitude  and  phase  variations  from  p'  to  p  on  path  m,  if  we 
suppress  the  time  delay  the  Extended  Huy  gens- Fresnel  Principle  becomes 


E'm(p,t)  -  dp'  Em{p' .t)ex^p')+i^p^ 


r\  he 


U-o(I+!p-p'|2/2L) 


HhxL 


(6.2) 


We  can  usually  take  \m  and  <pm  to  be  jointly  Gaussian,  allowing  us  to  construct  a  structure 
function7  [20] 


Dm{p-  p')  =  {(.\m(P0  +  P  Pq  +  p')  -  Xm{pO ,  Po))2) 

+  (  ( 0  m  (  PO  +  P j  Po  +  P  )  —  (pin{P( ):  Pq))~)  (6-3) 


which  satisfies  [17] 
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(6.4) 


for  the  refractive  index  structure  constant  C^rn(z)  =  10 ~12Cj-(z).  When  considering  prop¬ 
agation  that  is  approximately  parallel  to  the  earth,  it  is  reasonable  to  assume  that  the 
structure  constant  is  stationary,  i.e.  C.~  rn(z )  — >  C%>m.  Since  Xm  and  <pm  are  jointly  Gaussian 
we  can  define  a  complex  random  process 


V'm(P.P')  =  Xm(p-P')  +  i<t>m(p,p') 


(6.5) 


7The  Gaussian  assumption  and  its  associated  structure  function  are  limited  in  validity  to  what  is  known  as 
the  weak- perturbation  regime.  However,  a  more  general  derivation  will  lend  greater  validity  to  the  correlation 
function  expression,  given  below  in  Eq  (6.6),  that  will  suffice  for  our  purposes  in  this  chapter. 
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whose  correlation  function  on  path  m  can  be  evaluated  as 


-  exp  {-D(pi  -  Pi,  p2  -  pA)/ 2) 

Pi  ~  Pi?  +  (Pi  -  Pi )  ‘  (p2  -  Pi)  +  |p2  ~  Pi\2\  ,n 

% - )  (M>) 

where  we  have  used  the  square-law  approximation  to  the  rigorous  5/3-law  behavior  and 

Pm  -  {imh-ClmL)-3/s  for  m  =  R,  S,  T,  (6.7) 

is  the  turbulence  coherence  length  on  path  m,  under  the  assumption  of  a  constant  (7/  m 
profile.  In  what  follows,  we  shall  assume  that  there  is  sufficient  physical  separation  between 
the  reference,  signal,  and  target-return  paths  so  that  their  turbulence  effects  are  uncorrelated. 


6.2  Classical  Source 

6.2.1  Resolution  analysis 

We  start  by  looking  at  the  ghost  imaging  system  that  we  developed  in  Chapter  2,  which  uti¬ 
lized  a  pseudothermal  source  and  a  50  50  beam  splitter  to  create  signal  and  reference  fields 
with  a  phase-insensitensive  cross  correlation;  we  now  add  statistically  uncorrelated  turbu¬ 
lence  on  all  three  paths.  We  look  at  an  ensemble  average  of  the  correlation  of  photocurrents 
produced  by  each  detector,  using  the  same  assumptions  employed  in  our  resolution  analysis 
in  Chapter  2  to  arrive  at 

<cVi»  =  <ii(i)i2(f)> 

=  q2V2A i  JdTi  JdT2h(t-Tl)h(t-T2)  J  dp' t,)EJ (p',  T2)Egp',  r2)). 

Using  Eqs  (6.2)  and  (6.5)  to  back-propagate  /•.'•)  (//.  C)  to  E's(p2.t2)-  once  again  suppressing 
the  time  delays,  and  recognizing  the  independence  of  the  randomness  in  the  target  surface, 
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field  speckle,  and  turbulence,  we  have 


(E'Z{p1,T1)E'R(p1,T1)E2(p',  t2)E't(P\t2))  = 


k2 


4tt2  L2 


IdP2I 


dp:>>  e 


-ik0(\p' -Pt\2  +  \p' -p3f)/2L 


x  (e»K->'«)e«»'-»'){£;“(p1,Tl)£;yp1,r1)^(p2,T-2)Bi(p3,r2)>(r(p2)r(p3)>. 


Using  Eq  (2.5)  to  evaluate  the  surface  moment,  this  simplifies  to  to 


(i3i>,,n)£f(p,.r,)fi(p1,il)£'T(p'.T,)) 

=  p  Jip2  T(p2)(E%(pR  n)E's’(p2,  t2)E'r(Pi  .  Ti)  E'slp2-  t2)) 

Interestingly,  since  _  \  we  find  that  the  turbulence  on  the  target-return 

path  does  not  affect  the  resolution  of  the  final  image. 


In  our  resolution  analysis  in  Chapter  2  we  could  stop  here  and  use;  the  Gaussian-Schcll 
model  to  evaluate  the  fourth-order  field  moment.  However,  since  we  arc  allowing  for  turbu¬ 
lence  on  the  signal  and  reference  paths,  we  have  to  back  propagate  these  fields  to  the  source, 
giving  us 


{Er.{P1'  e)E's  (P2-  <2)E'R(pi,  Ti)E's(p2,  To))  —  J  ^Pi  Jdpo  J dPi  j dp. 

,  krP~iko{-L+\pl~P\\2l2L) 

<FB(PU)^>V  - 


l-  np-iko{L+\p2-p'2\2 /2L) 

EM.  _i! - - - 

42-kL  S  '  —i2nL 

„  hnPiko{L+\Pi-p'{\2 l2L) 

x  ER(p".  - Es(p't,  r2)e <*«■ 


J  ?  (  7”,— 1—  I  />«■»  - 


We  can  now  use  Gaussian  moment  factoring  and  Eq  (2.6)  to  evaluate  the  fourth-order  field 
moment,  and  Eq  (6.6)  to  evaluate  the  turbulence.  Taking  into  account  our  DC-block  filter 
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and  using  the  far-held  assumption,  we  have 


(Er(Pi;Ti)Es(p2,t2)E'r(p1,  T\)E's{p2,  t2))  = 


h  J  dpi  J  dpi  J  dpi 


x  e^pv[-p'\-p”) e^-P2-{P2-p'i) e-\p\~p”\2 EpRe-\p'2-p'i? /24 e 
x  e~\p\~P2\2 l2Poe-\p'2-P\\2 /2pI e-(T-2-n)2 /T$ 


(i^h+i^p+iPi'h+Ki2)/^ 


(6.8) 


In  evaluating  Eq  (6.8)  we  will  assume  that  the  turbulence  coherence  length  on  each  path 
m.  will  stay  large  enough  that  prn  »  p0.  as  p0  is  typically  on  the  order  of  a  few  wavelengths. 
We  will  still  assume  that  a0  p0,  but  we  will  make  no  assumption  about  the  relative  sizes 
of  «o  and  pm.  Under  these  conditions,  we  hnd  that 


mpi)) 


q2rfA x/f2  /  2 P 


2  pIpI 


Ipi 


2  2 


L 2  V  /  2/4/4  +  al  (pR  +  Ps) 


dp2T(p2)e  pl  mps+a>R+ps 


(6.9) 


Our  PSF  has  now  widened,  with  the  resolution  transforming  as  pL  -4  p'L  for 


Pl 


Pip 


2 PrPs  +  ao{p[ 


R 


Pl) 


2 PrPs 


(6.10) 


By  means  of  Eq  (6.7)  p'L  can  be  rewritten  purely  in  terms  of  source  size  and  structure 
constants  on  each  path  as 


Pl  = 


2  L 
k(l0 


1  +  22(1.09FL)«/5  (C^fs  +  (ClR) 


.6/5 


(6.11) 


It  is  also  worth  noting  that  turbulence  on  the  signal  and  reference  paths  have  identical 
impacts  on  the  PSF,  and  thus  identical  contributions  to  resolution  degradation. 


54 


6.2.2  Special  cases 

We  would  now  like  to  look  at  a  few  special  cases  that  illustrate  most  of  the  turbulence 
behavior  that  we  are  interested  in.  These  arc:  no  turbulence  on  either  path:  turbulence 
on  only  one  path,  which  corresponds  to  the  computational  case:  and  symmetric  turbulence 
on  both  paths,  which  is  a  good  approximation  when  both  paths  are  going  through  similar 
atmosphere. 


No  Turbulence  As  the  turbulence  coherence  lengths  increase  without  bound  (pR  -4 
oo.  ps  -4  oo)  the  effects  of  turbulence  vanish.  In  this  limit  p'L  =  pi-  and  Eq  (6.9)  be¬ 
comes 


(C(pi)) 


—Adi  A 


TV 


V- 


2  P 


7T  &  r 


dp2  Tip-)  ] 


.  i  cr/4 


(6.12) 


which  matches  Eq  (2.15).  our  result  calculated  without  turbulence. 


Computational  Case  When  doing  the  computational  case  the  reference  arm  is  calculated, 
so  there  is  no  turbulence  to  account  for  on  the  reference  path.  Thus  we  can  let  pR  -4  oo. 
and  we  find  that 


<<?(Pi)>  = 


2  2 
Q  T 


A,  A,  (2 P 


ncr, 


2p| 


2  Pi 


dp2  T(p 


•2)e 


\P\~P2\  2f)S 

4  a4+a  o 


(6.13) 


Be  comparing  Eqs  (6.9)  and  (6.13),  we  see  that  computational  and  non-computational  ver¬ 
sions  have  different  PSFs  and  that  the  resolution  for  the  computational  case  is  better.  That 
is, 


2p|  +  flu  „  I2PrPs  +  ao (Pr  +  Ps) 

2p|  PtV  ~ . “ 


(6.14) 


It  should  be  noted  that  while  this  situation  describes  the  computational  case,  it  is  possible 
to  have  a  non-computational  case  in  which  there  is  only  turbulence  on  one  arm.  Since  the 
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turbulence  on  both  paths  have  identical  effects,  a  situation  where  there  is  turbulence  on 
only  the  reference  path,  and  not  the  signal  path,  would  be  described  by  Eq  (6.13)  with  pR 
substituted  for  ps- 


Symmetric  Case  In  many  situations  the  signal  and  reference  paths  will  be  going  through 
atmosphere  with  similar  turbulence,  and  thus  approximately  the  same  structure  functions 
and  coherence  lengths.  This  lets  us  say  pR  =  ps  =  pu  and 


<C(p,»  =  ^d) —  f-^)  -rqH 

L-  \naLJ  Pt+OoJ 


\P\~P2\ 


(6.15) 


6.3  Non-classical  Source 


Having  looked  at  the  effects  of  turbulence  for  a  classical  light  source,  we  once  again  turn  our 
attention  to  non-classical  sources.  We  again  perform  a  Gaussian-state  analysis,  following  the 
resolution  derivation  in  Chapter  3  with  turbulence  added  on  all  three  paths.  For  propagation 
through  turbulence  we  will  replace  Eq  (6.2)  with  its  operator  equivalent 


p-p'\2/2L) 


i2nL 


for  to  =  S,  R.  T, 


where  ET(p.t )  —  E's(p,t)T{p)  +  Evac_s{p,  t)y/ 1  —  |T(p)|2  and  js  the  same  com¬ 

plex  random  process  from  Eq  (6.5)  which  encapsulates  the  effects  of  turbulence  on  path  m 
from  f)  to  p.  We  take  the  ensemble  average  of  the  photocurrents,  which  after  normally 
ordering  the  held  operators  and  applying  Eq  (3.6)  becomes 


(C(Pl))  =  (h(t)l2(t)) 

fdn  fdT2h(t-Tl)h(t-T2)  [  dp'(i^(p1:/ri)£r(A>,172)^(pli'r1)£^.(p',T 2)). 
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We  back  propagate,  evaluate  the  surface  moment  with  Eq  (2.5)  and  the  turbulence  with  Eq 
(6.6)  to  get 


(E'i(p  1,  Tl)E}(pf,  T2)E'R(Pl,  T\)E'T(p' .  r2)> 

=  J2  f  iP2T(p^){,El'R(pi.Tl)Ens(p2,T2)b'R(pi.T{)&s(p2,T2)). 


(6.16) 


Again,  we  find  that  there  is  no  effect  from  the  turbulence  on  the  target-return  path.  Propa¬ 
gating  back  to  the  source,  evaluating  the  field  and  turbulence  moments,  and  implementing 
the  DC-block  filter,  we  arrive  at 


(Er{PU  Tl)E^(p2,  To)  E,R{p1,T1)E's(p2,  T2))  = 
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(6.17) 


Equation  (6.17)  can  he  directly  evaluated,  and  we  find  the  photnenreent  average  to  he 
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(6.18) 


where  X  =  XT0Po/ao  —  XToPl/a|  is  the  brightness  term  from  the  SNR  derivations.  Since 
wc  are  operating  in  the  far  field,  there  is  no  resolution  gain  from  the  use  of  entangled  signal 
and  reference  fields,  even  in  the  low-brightness  limit  in  which  the  Gaussian  state  becomes 
the  biphoton  state.  Turbulence  causes  the  same  resolution  spreading  in  phase-sensitive 
ghost  imaging  as  it  does  in  phase-insensitive  ghost  imaging,  regardless  of  the  nature  of  the 
correlation  (classical  or  non-classical). 
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6.4  Laser  Radar 


Finally,  we  would  like  to  compare  the  turbulence-induced  degradation  suffered  by  ghost 
imaging  systems  to  that  for  a  laser  radar  system  under  the  same  atmospheric  conditions. 
We  will  be  using  the  laser  radar  system  developed  in  Chapter  5.  and  start  by  taking  an 
ensemble  average  of  the  image-bearing  term  in  Eq  (5.1)  to  arrive  at 


(■ I(Pi ))  =  qvAi{EZ(p!,t)E3(pi,t)). 


(6.19) 


Propagating  back  to  the  lens,  we  have 


(I(pi))  =  qrjAi  J  dp'  J  dp"(E2{p',t)E2{p",t)) 

A'j  Ai 
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—i2ixL 


i‘2nL 


(6.20) 


We  are  not  considering  any  turbulence  between  the  lens  and  the  CCD  array.  For  con¬ 
venience  we  are  using  a  1:1  imaging  system,  but  in  a  realistic  implementation  the  image 
produced  is  a  minified  version  of  the  target.  For  this  system  we  are  considering  turbulence 
on  both  paths,  viz.,  on  the  signal  path  (S),  from  the  source  to  the  target,  and  the  target- 
return  path  ( T ).  from  the  target  to  the  lens.  Thus,  back  propagating  to  the  target  we 
have 


(E*(p':t)E2(p":t))  ^Jdp2J dp3(e^^^^'’^)(T*(P2)T(p3)) 
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Using  Eq  (2.5)  and  the  far-held  assumption  this  simplifies  to 
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Now.  back-propagating  (El(p2,t)Ei(p2.t)),  and  assuming  the  far-held  condition,  we  have 
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Our  laser  radar  system  is  constrained,  for  comparison  purposes,  to  have  an  on-target  average 
illumination  pattern  equal  to  that  of  the  ghost  imager.  Since  the  laser  is  a  coherent  source, 
that  means  that  the  source  beam  waist  wQ  for  this  system  is  approximately  the  same  size  as 
the  source  coherence  length  p0  for  the  ghost  imaging  system.  As  such,  we  can  assume  that 
for  the  laser  radar  system  wq  <  ps,  which  let  us  simplify  to 

.9 

2p  2\P2\ 

(El{p2.t)Ei(p2,  t))  —  -e  "l  , 

Txap 


indicating  that  for  this  system  there  is  no  effect  on  the  resolution  from  turbulence  on  the 
signal  path. 


This  lets  us  rewrite  Eq  (6.20)  as 


A'2  A-2 


-|p"-p'|2/2p2e-Ta(p"-p')'Pii 


(6.21) 


where  T(£)  is  again  the  two-dimensional  spatial  Fourier  transform  of  T{p)-  We  invert  the 
coordinates  as  pi  — >  —pi  so  that  we  have  an  upright  image,  and  use  difference  coordinates 
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to  evaluate  the  integrals  over  the  lenses  to  arrive  at 
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To  evaluate  the  spreading  due  to  turbulence  we  can  approximate  the  Jine  PSF.  which  has  its 
first  zero  at  \px\  —  1.22A 0L/D2i  as  a  Gaussian  PSF  with  an  e_1  point  at  \px  \  =  2XqL/hD2 , 
i.e.. 
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This  lets  us  simplify  Eq  (6.22)  to 
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We  can  now  say  that  the  effect  of  turbulence  on  the  resolution  is  to  degrade  it  by  a  factor 
of  \j “^2 pi,  For  a  comparison  to  the  ghost  imaging  system,  we  once  again  say  that  all 
lens  diameters  are  of  the  same  size,  so  that  D  =  2%,  and  our  spreading  factor  becomes 
yj  2pi  Q°  •  This  is  the  same  spreading  we  found  when  there  was  turbulence  on  one  arm  of  the 
ghost  imaging  system,  as  shown  in  Eq  (6.13).  From  this  wo  can  conclude  that  ghost  imaging 
systems  can  Ire  constructed  that  do  not  suffer  resolution  degradation  from  turbulence  worse 
than  the  degradation  found  in  a  flood-light  illumination  laser  radar  system,  ft  should  also  be 
noted  that  for  a  raster-scanning  type  laser  radar  system,  which  we  have  not  yet  considered, 
the  source  beam  size  is  significantly  larger  than  for  the  flood-light,  system,  and  therefore 
the  turbulence  on  the  signal  path  will  no  longer  be  insignificant.  Indeed,  partially  coherent 
beams,  such  as  our  pseudothermal  source,  have  been  shown  to  be  less  affected  by  turbulence 
than  similar  fully  coherent  beams  [21,  22],  indicating  that  ghost  imaging  might  have  an 
advantage  in  comparison  to  a  raster-scanning  laser  radar.  We  also  might  expect  that  the 
effects  of  turbulence  on  the  SNR  of  ghost  imaging  and  laser  radar  systems  could  be  different, 
but  we  have  not  explored  this  issue. 
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Chapter  7 


Improving  Ghost  Image  Quality 


So  far  we  have  analyzed  reflective  ghost  imaging  for  classical  and  non-classical  sources,  as 
well  as  for  a  computed  reference  beam.  These  results  were  compared  to  a  laser  radar  system; 
ghost  imaging  has  a  somewhat  advantageous  asymptotic  SNR  behavior  for  high-brightness 
illumination  and  long  capture  times,  as  well  as  improved  low-brightness  behavior  when  there 
is  a  non-classical  correlation  between  the  fields,  and  when  a  computed  reference  arm  is 
employed.  However,  the  SNR  is  still  limited  by  the  speckled  nature  of  the  return  patterns 
cast  by  the  rough  surface  of  the  target.  In  this  chapter  we  propose  three  methods  to  improve 
the  image  quality,  focusing  on  reducing  the  target-speckle  imposed  SNR  limit. 

First,  we  propose  a  detection  scheme  in  which  multiple  bucket  detectors  are  employed, 
each  capturing  the  returned  light  at  different  transverse  locations  within  the  same  plane,  as 
detailed  in  Fig.  7-1.  The  target-speckle  attributes  of  the  returning  held  decorrelate  with 
transverse  location,  so  each  detector  should  see  different  speckle  behavior,  allowing  us  to 
average  it  out  and  improve  the  asymptotic  SNR.  This  same  target-speckle  averaging  could 
also  be  accomplished  by  the  use  of  multiple  wavelength  illumination,  as  shown  in  Fig.  7-4.  If 
we  use  several  narrow-band  sources  sufficiently  far  apart  in  wavelength,  the  return  patterns 
from  each  source  will  have  uncorrelated  speckle  statistics  at  the  same  transverse  location, 
allowing  us  average  out  the  speckle  effects.  Finally,  each  bucket  detector  measurement  in 
a  ghost  imaging  setup  can  be  thought  of  as  a  random  projection  of  the  target  reflection 
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pattern.  Since  we  know  the  pattern  being  cast  on  the  target,  compressive  sensing  methods 
can  be  applied  to  the  bucket  detector  measurements  to  achieve  a  high-quality  image  quickly. 
This  connection  has  already  been  made,  and  experimentally  verified,  in  [12]  for  transmissive 
ghost  imaging.  In  reflective-imaging  we  have  the  added  complication  of  the  target-speckle, 
but  we  will  demonstrate  that  a  connection  can  still  be  made  back  to  the  idea  of  projective 
measurements,  enabling  the  use  of  compressive  sensing  techniques. 


7.1  Multiple  Bucket  Detectors 


For  this  setup  we  consider  n  co-plan ar  bucket  detectors  in  some  arbitrary  setup.  The  corre¬ 
lation  function  at  each  transverse  location  px  is  taken  to  be  the  average  of  the  correlation 
function  for  the  CCD  pixel  at  p\  and  each  bucket  detector.  The  output  photocurrent  of  a 
CCD  pixel  will  still  be  denoted  as  ii(t).  while  the  output  of  each  bucket  detector  i  will  be 
labeled  i({t)  for  2  <  f:  <  n  +  1.  We  thus  build  up  the  correlation  function  pixel- wise  as 
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C(P l)  =  -  XI 


i  rT' 12 


£—2  I  **  “2  /  /2 


dt  i\  (t)ic(t). 


(7.1) 


The  resolution  and  contrast  for  this  system  can  easily  by  found  from 


1  n+1  i  fTj/2 

(C(p i))  =  -£t  / 

n  /=9  J —7'ji 


'll  2 


(7.2) 


by  recognizing  that  (ii(t)ie(t))  is  the  same  for  all  £,  and  not  a  function  of  time.  The  average 
correlation  function  becomes  (C(pi))  =  (ii(t)ie(t)),  which  is  what  was  derived  in  Chapter  2 
for  one  bucket  detector.  Thus,  no  resolution  or  contrast  gain  is  realized  by  utilizing  multiple 
bucket  detectors. 

However,  this  setup  does  offer  a  SNR  enhancement.  Using  the  SNR  definition  in  Eq 
(2.21),  as  well  as  the  same  simplifying  assumptions  that  allow  us  to  use  Eq  (2.22)  for  the 
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Bucket  -  ►  Correlator  ◄ 


Figure  7-1:  Reflective  ghost  imaging  with  multiple  bucket  detectors. 


first  moment,  we  are  left  only  needing  to  calculate  the  second  moment  of 


(C2(p  i)) 


n  +  Wi+l  .  rT,/2  rT,/  2 


(7.3) 


We  are  interested  in  raising  the  high-brightness,  long  integration  time  SNR  asymptote,  so 
we  restrict  our  analysis  to  this  limit  ,  in  which  we  find 


SNR„.  long-T; 


n22A2l'Ka2 


n+1  n+1 

EEr< 


C= 2  m=2 


(7.4) 


where 


J  du  e-|l/-A<m|2/20(i/.  4a). 


(7.5) 
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The  T  function  in  Eq  (2.26)  conies  from  integrating  over  the  bucket,  detector  in  two  different 
coordinates  systems  p'  and  p"  a  function  of  the  difference  of  these  coordinates.  To  evalu¬ 
ate  it  we  moved  to  the  normalized  difference  coordinates  v  =  pzkoip'  —  p")/L  and  integrated 
over  the  two-circle  overlap  function.  This  worked  because  the  different  coordinate  systems 
were  for  the  same  bucket  detector,  and  centered  at  the  same  point.  Since  we  are  now  using 
multiple  bucket  detectors,  we  have  to  account  for  the  difference  in  their  transverse  location. 
We  want  to  recenter  one  coordinate  system  so  that  both  are  aligned  by  shifting  it  by  some 
p[  —  p^  where  p'(  is  the  center  of  detector  £.  and  p"  is  the  center  of  detector  m.  Thus,  we 
now  define  v  =  pLk0(p'  -  p") / L  +Ac,m  for  Atm  =  pLk0(p'e-  p"n)/L,  allowing  us  to  correctly 
use  the  two-circle  overlap  function  to  evaluate  our  expression. 


The  averaging  over  speckle  statistics  is  encapsulated  in  T t.m,  and  is  where  we  sec  the  ad¬ 
vantage  of  using  multiple  bucket  detectors.  As  the  distance  between  the  detectors  increases. 
I  \.m  and  the  deleterious  speckle  effects  decrease,  which  leads  to  an  SNR  increase.  Looking 
again  at.  detector  sizes  of  A2  —  Trajp  and  placing  two  detectors  side-by  side  (centers  separated 
by  2a0).  we  find  T t.m  —  0.049.  which  is  a  significant  decrease  from  I'c.r  =  T  =  .61:  this  con¬ 
figuration  yields  an  SNR  of  6.03.  If  we  separate  the  detectors  by  a  distance  of  2 a0  (so  that 
the  centers  are  4a0  apart)  we  have  T £.rn  —  4.3  x  10-7,  and  a  SNR  of  6.52.  which  is  twice  that 
of  the  single  detector  case.  Thus,  detectors  of  this  size  that  are  side-by-side  have  slightly 
correlated  speckle  statistics,  while  detectors  with  centers  separated  by  4a0  have  essentially 
uncorrelated  speckle  behavior,  as  can  be  seen  in  Fig.  7-2.  This  makes  intuitive  sense  because 
the  speckle  fluctuation’s  average  size  are  inversely  proportional  to  the  coherence  length  of 
the  light  illuminating  the  target  as  ‘speckle  length’  =  2 L/k0pL  =  a0.  Locations  separated 
by  lengths  significantly  greater  than  a0  in  the  bucket  detector  plane  will  have  uncorrelated 
speckle  statistics.  This  means  that  adding  bucket  detectors  far  enough  apart  linearly  in¬ 
creases  the  SNR.  while  detectors  placed  closer  together  will  yield  a  slightly  lower  increase, 
as  seen  in  Fig.  7-3.  This  speckle  averaging  is  the  same  effect  that  is  seen  when  increasing 
the  size  of  the  bucket  detector. 
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(a)  IYm  as  a  function  of  the  distance  between  (b)  SNR  as  a  function  of  the  distance  between 
the  center  of  detectors  (  and  in.  the  center  of  detectors  t  and  m. 

Figure  7-2:  VLm  for  distances  d  =  0  (t  =  m )  to  d  =  4a0  between  the  centers  of  detect  ors 
d  and  m,  and  the  subsequent  SNR  for  the  two  detector  configuration,  are  plotted.  The 
detectors  are  taken  to  have  area  A2  =  wa#. 


Figure  7-3:  The  SNR  per  bucket  detector  for  square  detector  arrays.  The  collecting  lens  for 
each  detector  is  taken  to  be  a  circle  of  area  A2  =  na^.  and  they  are  arranged  is  squares  (with 
sides  touching)  of  size  1.  2  x  2.  3  x  3,  4  x  4,  5  x  5 
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7.2  Multiple  Wavelengths 


In  place  of — or  possibly  in  conjunction  with — using  multiple  bucket  detectors  we  can  use  a 
multiple  wavelength  source  to  increase  the  SNR.  Whereas  the  bucket  detectors  in  Sec.  7.1 
see  different  speckle  statistics  due  to  their  transverse  separation,  the  return  fields  for  the 
various  wavelengths  have  different  speckle  statistics  at  the  same  transverse  location  because 
each  wavelength  diffracts  differently  from  the  rough  surface.  This  will  allow  us  to  make 
multiple  co-linear  measurements  to  average  out  the  effects  of  the  target  speckle. 

For  the  multiple  bucket  detector  setup  we  did  a  full  treatment,  finding  the  correlation 
between  the  speckle  statistics  for  bucket  detectors  separated  by  arbitrary  distances.  This 
was  because  it  is  reasonable  to  assume  that  in  an  implementation  we  could  be  constrained 
by  the  size  of  the  array:  we  might  want  to  pack  the  detectors  tightly.  We  should  not  be 
constrained  to  choosing  similar  wavelengths  for  our  sources.  Thus,  we  are  only  interested 
in  finding  the  cut-off  wavelength  separation  at  which  the  return  fields  have  uncorrelated 
target-speckle  statistics,  and  the  image  characteristics  in  this  regime. 

We  start  with  the  setup  in  Fig.  7-4,  and  we  build  up  the  correlation  function  as  the 
average  of  separate  correlation  measurements  made  at  every  wavelength.  Each  wavelength 
diffracts  differently,  so  if  we  assume  the  same  source  coherence  and  intensity  radii,  the  on- 
target  average  illumination  pattern  will  be  more  spread  out  for  some  wavelengths,  less  for 
others.  This  will  result  in  less  power  being  measured  at  each  detector  (both  CCD  and  bucket) 
at  certain  wavelengths;  therefore,  when  averaging  the  separate  correlation  functions,  we  need 
to  scale  the  measurements  appropriately.  For  our  work,  we  will  scale  each  measurement  by 
(Pl/ 444  V^n)2!  where  4°  and  pf*  are  the  on-target  intensity  and  coherence  radii  at 
wavelength  \e.  and  dL  and  pL  are  the  average  of  (44  and  {44  respectively.  Following 
our  previous  work,  CCD  output  44)  associated  with  wavelength  \t  is  correlated  with  the 
bucket  detector  output  lip  (t)  for  each  £  and  averaged  as 
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The  resolution  and  image  contrast  are  found  from  an  an  ensemble  average  as 


(C(P  i)> 


n  ' 


6=1 


rT,/2 


-27/2 


dt(i}{\t)iflt)). 


(7.7) 


Since  each  is  a  function  of  the  wavelength  Xe,  they  each  have  a  different  point- 

spread  function,  and  therefore  the  resolution  is  limited  by  the  longest  wavelength  used  to 
ma Xf(p^)  =  2L/\mni(ke)ao.  Looking  at  the  contrast,  we  find  the  background  term  becomes 
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q2rfAiA2 
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j dp2  T(p2); 
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and  the  image-bearing  term,  under  the  assumption  that  each  wavelength  is  short  enough  to 
resolve  the  target,  becomes 
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q2ifA1A2 
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(7.9) 


Using  the  definition  in  Eq  (2.16).  we  find 
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At 
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(7.10) 


This  differs  from  the  results  for  using  multiple  bucket  detectors  only  slightly,  as  the  coherence 
lengths  are  added  in  parallel:  since  the  oil-target,  coherence  lengths  the  will  be  of  the  same 
order  of  magnitude,  the  contrast  has  the  same  basic  interpretation  as  being  approximately 
the  inverse  of  the  number  of  resolution  cells  in  the  image. 


Turning  to  the  SNR  evaluation,  and  the  definition  in  Eq  (2.21),  we  are  left  only  needing 


to  evaluate  the  second  moment  of  the  correlation  function  as 
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(7.11) 


We  again  look  at  the  high-brightness  long-integration-time  limit  where  the  effects  of  the 
target  speckle  dominates.  To  do  this  we  need  to  review  our  treatment  of  the  scattering 
effects  ol  the  rough  surface  given  in  Eq  (2.4).  This  treatment  was  for  a  single  wavelength 
A0.  and  the  phase-insensitive  autocorrelation  is  given  by  Eq  (2.5).  We  are  now  dealing  with 
multiple  wavelengths  of  light,  so  we  define  the  target  reflection  coefficient  for  wavelength  X( 

cLS 


Te(p)  =  AnP)e2lklA:ip):  (7.12) 

where  Ay  —  2tt/A£.  The  phase-insensitive  cross  correlation  of  reflection  coefficients  for  differ¬ 
ent  wavelengths  is  thus 

C (p'))  -  v^Vj(e'2i(^;(phfemA;(p))) 

—  sJAATAp' )e"  2<T?(fc?+b27  ,)+^ki:k"iKzz{p-p') 

As  before,  we  are  assuming  that  the  covariance  function  Kzz(p  —  p')  for  the  height  variations 
Az(p)  is  very  narrow,  allowing  us  to  approximate 

C Tiip)Tm(p '))  =  A, A mS(p  -  mvxr*9*-*’#  (7.13) 


Thus,  when  n:(%  -  km)2  =  1,  then  (T(*(p)Tm(p'))  =  A,AmA(p  -  p')T(p)e~-  «  0.  and 
the  return  fields  from  the  two  different-wavelength  sources  are  essentially  uncorrelated.  This 
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Figure  7-4:  Reflective  GI  system  using  multiple  lasers  at  different  wavelengths,  and  one 
bucket  detector  for  each  wavelength. 


gives  us  the  minimum  wavelength-separation  criterion  between  two  wavelengths  as 


AA|  > 


A2 

27TCT, 


(7.14) 


where  A  is  approximately  the  average  of  the  wavelengths.  In  this  regime,  we  can  approximate 
the  exponential  as  a  Kronecker  delta  function,  i.e.,  e--az(ki-kmY  5(m.  This  gives  us  a  high¬ 
brightness  long-integration-time  SNR  of 


SNR//  long-2'/ 


n2,42 

7raoF 


(7.15) 


This  is  the  same  basic  result  we  found  for  the  multiple  bucket  detector  setup,  when  the 
bucket  detectors  were  sufficiently  separated;  averaging  over  multiple  correlations,  each  with 
uncorrelated  speckle  statistics,  leads  to  a  linear  increase  in  the  maximum  SNR. 


7.3  Compressive  Sensing 


In  traditional  sensing  methods,  resolution  is  limited  by  the  Nyquist  rate,  i.e. ,  you  must 
sample  at  a  rate  at  least  twice  as  high  as  the  desired  signal  bandwidth  in  order  to  acquire 
them.  However,  when  dealing  with  signals  that  have  a  sparse  representation  in  some  basis, 
the  captured  data  can  then  be  compressed,  reducing  the  number  of  bits  needed  to  accurately 
represent  it  [13].  A  sparse  representation  means  that  in  some  basis,  almost  all  coefficients 
arc  close  to  zero,  so  most  of  the  information  is  held  in  only  a  few  significant  coefficients  [14]. 
For  a  simple  example,  imagine  some  audio  signal  composed  of  only  three  frequencies.  If  we 
wanted  to  sample  in  time,  then,  according  to  Nyquist,  we  would  have  to  sample  at  twice  the 
highest  frequency  to  have  a  faithful  representation  of  the  signal.  However,  in  the  frequency 
domain,  this  signal  can  be  represented  by  only  three  coefficients;  if  we  wanted  to  compress 
the  signal,  we  could  represent  it  in  the  frequency  domain  and  only  keep  the  three  important 
data  points  corresponding  to  the  three  frequencies  present,  significantly  reducing  the  number 
of  data  points  required  to  represent  the  signal. 

Any  signal  with  some  definable  structure  is  likely  to  have  a  sparse  representation  in 
some  basis  [14].  Of  interest  to  us  are  "natural  images."  Natural  images  are  generally 
discontinuously  smooth;  that  is.  they  are  composed  of  areas  of  slowly  varying  features  that 
are  separated  by  sharp  boundaries.  As  is  exploited  in  a  number  of  image  compression 
schemes,  such  as  the  ubiquitous  JPEG  format,  natural  images  have  a  sparse  representation 
in  the  2-D  discrete  cosine  basis  [14], 

Compressive  sensing  is  a  means  to  directly  acquire  the  significant  coefficients  of  a  signal 
in  its  sparse  basis,  requiring  significantly  fewer  measurements  than  if  the  signal  was  sampled 
at  the  Nyquist  rate  [13].  Conceptually,  this  technique  exploits  our  knowledge  of  the  structure 
of  the  image;  we  know  that  it  is  sparse  in  some  basis  (even  though  we  don’t  know  which 
coefficients  matter),  and  use  this  knowledge  to  reduce  the  amount  of  information  we  need  to 
acquire. 

At  the  heart  of  this  technique  is  the  notion  that  the  information  is  preserved  through 
linear  transformations:  therefore,  by  taking  a  series  of  linear  projections,  enough  information 
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can  be  gathered  to  reconstruct  the  signal.  This  process  involves  projecting  the  sparse  target 
onto  a  series  of  vectors  to  build  up  a  vector  of  projections  [14].  This  series  of  projections 
takes  the  form  of  a  measurement  matrix  A  such  as 

y  =  Ax  (7.16) 

where  a:  is  a  sparse  representation  of  our  signal,  and  y  is  the  vector  of  projections.  The 
most  efficient  measurement  matrices  have  the  lowest  correlation  between  columns  of  the 
measurement  matrix.  This  is  often  referred  to  as  incoherence  [14],  and  means  that  for  each 
projection,  each  coefficient  of  x  receives  a  nearly  uncorrelated  weighting.  To  recover  x  one 
uses  linear  programming  methods  and  knowledge  of  the  measurement  matrix  to  find  x  based 
on  an  Li-optimization. 

The  authors  of  [12]  realized  that  the  output  of  a  bucket  detector  in  transmissive  ghost 
imaging  is  simply  the  projection  of  the  target  transmissivity  pattern  onto  some  illumination 
pattern.  Therefore,  by  sending  a  series  of  illumination  patterns  at  the  target  a  data  set  of 
minimally  correlated  projections  can  lie  built  up.  which  is  the  ideal  data  set  for  compressive 
sensing.  By  using  an  SLM,  they  controlled  the  illumination  patterns  cast  on  the  target,  and 
used  that  knowledge  to  perform  the  reconstruction.  In  traditional  ghost  imaging,  the  known 
illumination  pattern  is  used  to  perform  a  correlation  with  the  bucket  detector’s  output  to 
form  an  image.  However,  by  instead  performing  compressive  sensing  on  that  output  (again 
using  knowledge  of  the  illumination  pattern)  the  sparsity  of  natural  images  is  leveraged  to 
perform  a  better  reconstruction  from  fewer  measurements. 

The  projective  measurement  made  by  the  bucket  detector  in  transmissive  ghost  imaging 
can  be  viewed  in  terms  out  the  output  current  as 

h(/)  '  [  dpl{p.t)\T(p)\\  (7.17) 

JAz 

where  I(p,  i)  =  \E(p,t)\2.  and  we  are  ignoring  shot-noise  for  the  moment.  This  can  be 
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written  in  matrix  form  as 


i>  =  IT  (7.18) 

where  T  is  a  column  vector  of  |T(p)|2  indexed  by  the  transverse  location  p.  I  is  the  pseudo¬ 
random  measurement  matrix  with  columns  indexed  by  time  t,  and  rows  indexed  by  p.  Our 
output  is  thus  a  column  vector  indexed  by  time.  Now,  T  is  not  sparse  in  its  spatial  repre¬ 
sentation,  but  for  the  rotation  matrix  to  some  sparse  basis,  we  have 

T  —  'E'T'  (7.19) 

where  T'  is  our  representation  of  the  target  transmissivity  in  the  sparse  basis.  If  we  combine 
the  basis  rotation  matrix  and  measurement  matrix  as  A  —  J\Er  we  then  have 

i2  =  AT'.  (7.20) 


In  reflective  imaging,  it  is  a  little  more  complicated  to  frame  the  bucket  detector's  output 
as  a  projective  measurement,  due  to  the  scattering  from  the  rough  surface.  However,  we  can 
make  the  connection  as  follows.  First,  since  the  bucket  detector  measures  the  total  power 
that  passes  through  the  collecting  lens,  we  can  theoretically  look  in  any  plane  behind  the 
lens  to  make  our  connection:  the  total  power  will  be  the  same  in  every  plane.  Therefore, 
we  can  look  in  the  image  plane  for  our  connection,  where  the  transverse  intensity  pattern 
is  a  product  of  the  illumination  pattern  I(p.t)  =  \E's(p,t)\2,  the  target  reflectivity  T(p)t 
and  and  exponentially  distributed  random  process  S(p).  Here  S(p)  accounts  for  the  effects 
of  the  rough-surface  scattering,  and  for  each  pixel  in  the  final  image  can  be  thought  of 
as  a  collection  of  independent,  identically-distributed  (IID)  exponential  random  variables. 
Therefore,  we  can  write  the  output  of  the  bucket  detector  as 

*2 (t)=  [  dp  I(p.  t)S(p)T(p).  (7.21) 

Ja, 
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This  can  be  put  in  matrix  form  as 


i2  =  1ST  (7.22) 

where  I  is  the  intensity  matrix  illuminating  the  target  as  before,  T  is  a  column  vector  of 
T  indexed  by  transverse  location  p,  S  is  a  diagonal  matrix  with  each  element  being  an  III) 
exponential  random  variable,  and  i2  is  our  projective  measurements,  indexed  by  time.  Thus, 
the  only  difference  between  the  transmissive  and  reflective  cases  is  the  diagonal  matrix  S. 
Now.  rotating  to  a  sparse  basis  as  in  Eq  (7.19),  have  an  output 

i2  =  AT1  (7.23) 


for  A  =  I  Sty.  While  our  measurement  matrix  still  has  columns  that  are  fairly  uneorrelatod, 
we  no  longer  know  the  measurement  matrix.  Whether  we  can  fully  recover  V ,  or  merely  a 
speckled  version,  is  the  subject  of  current  research. 
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Chapter  8 


Conclusion  and  Future  Work 


In  this  thesis  we  developed  a  framework  for  analyzing  reflective  ghost  imaging  systems  in  the 
context  of  stand-off  sensing.  We  developed  results  for  three  possible  imaging  systems,  which 
were  then  compared  to  each  other  as  well  as  to  a  basic  laser  radar  system.  We  then  applied 
turbulence  on  the  propagation  paths  for  these  systems  and  analyzed  the  resulting  resolution 
degradation.  Finally,  we  explored  ways  to  improve  reflective  ghost  imagers,  primarily  the 
high-brightness,  long  integration  time  SNR. 

The  first  system  we  explored  utilized  pseudothermal  light  for  the  signal  and  reference 
beams,  with  a  classical  phase-insensitive  cross  correlation  between  them.  This  is  realized  by 
passing  laser  light  through  rotating  ground  glass,  and  then  a  50  50  beam  splitter  to  create 
identical  spatiotemporally  random  fields  on  both  arms.  This  system  is  perhaps  the  most 
robust  of  those  wc  explored;  the  laser  source  can  be  made  very  strong  without  affecting  the 
correlation  between  the  fields,  and  wc  do  not  need  to  worry  about  the  calibration  of  an  SLM. 
The  next  system  we  investigated  utilized  the  output  of  an  SPDC,  along  with  a  polarizing 
beam  splitter,  to  create  entangled  signal  and  reference  beams  with  a  non-classical  phase- 
sensitive  cross  correlation.  We  evaluated  these  fields  within  a  Gaussian-state  framework, 
where  the  low-brightness  limit  is  the  biphoton  state  and  the  high-brightness  limit  is  classically 
correlated  light.  Finally,  we  investigated  a  computational  ghost  imaging  system  in  which  the 
signal  arm  is  created  by  modulating  the  phase  front  of  a  laser  with  an  SLM,  and  the  reference 
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field  is  calculated  from  the  known  applied  phase  modulations.  The  laser  radar  system  we 
used  for  comparison  is  a  flood-light  system  in  which  the  target  is  uniformly  illuminated,  and 
a  CCD  array  is  used  to  make  a  high-resolution  spatial  measurement  of  the  reflected  field. 

We  compared  these  systems  through  their  resolutions,  contrasts,  and  signal-to-noise  ra¬ 
tios.  To  facilitate  the  comparison  we  said  all  lenses  were  of  identical  area  ttoq,  where  a0  is 
the  average  intensity  radius  of  the  source  for  ghost  imagers.  We  found  that  all  three  ghost 
imagers,  as  well  as  the  laser  radar  system,  have  the  same  effective  resolution  of  pL  =  2L/k0a0. 
Looking  at  the  contrast,  the  pseudothermal  ghost  imager,  SPDC  ghost  imager  in  the  high¬ 
brightness  regime,  and  the  computational  ghost  imager  have  a  constant  background  term 
which  limits  the  contrast  to  Eq  (2.18).  This  can  be  overcome  by  AC-coupling  the  outputs  of 
either  detector.  In  the  low-brightness  regime  the  output  of  the  SPDC  becomes  the  biphoton 
state,  and  the  resulting  images  are  background  free.  This  is  also  true  for  the  laser  radar, 
and  thus  both  systems  have  an  unlimited  contrast  as  defined  in  Eq  (2.16).  For  AG-coupled 
operation,  all  ghost,  imager  systems  and  the  laser  radar  system  have  the  same  resolution  and 
contrast.  However,  a  performance  difference  can  be  found  in  the  SNR  of  each  system. 

The  maximum  SNR  for  each  system  is  achieved  when  the  source  is  strong  and  the  inte¬ 
gration  time  is  long.  This  averages  out  the  shot  noise  in  the  detectors,  and  for  ghost  imaging, 
the  speckled  nature  of  the  spatiotcmporally  varying  light  that  illuminates  the  target.  How¬ 
ever,  the  roughness  of  the  surface  creates  speckle  in  the  return  field  that  is  time-invariant; 
integrating  for  longer  will  not  remove  this  randomness.  These  target-speckle  variations  arc 
different  at  each  transverse  location,  and  are  the  limiting  factor  in  the  SNR  for  both  ghost 
imagers  and  laser  radars. 

However,  ghost  imagers  have  an  advantage.  The  size  of  the  target-induced  speckle  is 
inversely  proportional  to  the  coherence  length  of  the  light  at  the  target.  Thus,  the  average 
size  of  these  speckles  for  a  ghost  imager  is  o0  =  2 L/k^p^.  If  our  receiving  lens  is  of  the  same 
size  as,  or  larger  than,  the  transmitting  lens,  there  will  be  some  averaging  of  the  speckle  by 
the  bucket  detector.  For  laser  radar  the  target  is  illuminated  with  nearly  uniform  light,  so 
the  speckle  size  cast  on  the  receiving  lens  is  much  smaller  than  for  a  ghost  imager.  However, 
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for  ghost  imaging  the  light  is  captured  by  a  bucket  detector,  while  for  a  laser  radar  the  light 
is  focused  onto  a  CCD  array;  the  speckle  size  is  much  larger  than  a  CCD  pixel,  resulting  in 
almost  no  speckle  averaging.  Therefore,  while  the  laser  radar  system  is  limited  to  SNR  =  1, 
the  ghost  imager  can  achieve  an  SNR  in  the  low  single  digits.  As  discussed  in  Chapter  2,  for 
a  receiving  lens  of  area  A2  =  vra^.  the  maximum  SNR  for  the  ghost  imagers  is  3.266. 

For  low-brightness  situations,  we  saw  a  larger  difference  in  the  SNR  behavior  of  the 
various  systems.  For  the  pseudothermal  ghost  imager,  we  make  two  measurements  of  classical 
light,  and  thus  have  randomness  from  both  measurements.  For  the  SPDC  imager,  the  low- 
brightness  limit  is  approximately  the  biphoton  state;  while  we  still  have  randomness  from 
the  two  measurements,  it  is  significantly  lower  than  for  the  pseudothermal  imager.  However, 
the  computational  imager  has  a  better  low-brightness  SNR  asymptote  than  both  systems, 
as  it  only  has  randomness  from  one  measurement.  The  only  area  where  the  laser  radar 
has  better  SNR  characteristics  than  the  ghost  imagers  is  for  very  short  integration  times, 
and  very  high-brightness  illumination.  This  is  because  the  laser  radar  only  has  to  overcome 
shot  noise,  which  can  be  accomplished  with  either  a  long  integration  time  or  high-brightness 
illumination;  the  ghost  imagers  also  have  to  average  out  the  randomness  of  the  source,  which 
can  only  be  accomplished  with  a  sufficiently  long  integration  time. 

After  exploring  the  basic  image  characteristics  for  these  systems,  we  looked  at  the  effects 
of  turbulence  on  the  image  resolution.  We  found  that  ghost  imagers,  classical  or  quantum, 
have  the  same  resolution  degradation  from  turbulence.  There  was  less  degradation  for  com¬ 
putational  imaging,  as  there  is  no  turbulence  on  the  reference  path.  It  is  also  worth  noting 
that  turbulence  on  the  return  path,  from  the  target  to  the  bucket  detector,  has  no  effect 
on  resolution.  In  comparison,  the  laser  radar  system  had  the  same  degradation  as  the  com¬ 
putational  imager.  Thus,  ghost  imagers  and  laser  radar  systems  suffer  similar  resolution 
degradation  from  turbulence. 

Finally,  we  looked  at  ways  of  improving  the  SNR  performance  of  ghost  imaging.  We 
showed  that  bucket  detectors  spaced  far  enough  apart  measure  light  with  uncorrelated  target- 
speckle  statistics:  thus,  by  using  multiple  bucket  detectors  at  different  transverse  locations, 
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we  can  increase  the  maximum  SNR  linearly  with  the  number  of  detectors.  We  also  showed 
that  light  sources  far  enough  apart  in  wavelength  result  in  uncorrelated  target-speckle  statis¬ 
tic  at  the  bucket  detector:  thus,  by  measuring  each  wavelength  separately  we  can  increase 
the  SNR  linearly  with  the  number  of  source  wavelengths  utilized.  We  also  explored  applying 
compressive  sensing  methods  to  improve  image  quality  while  reducing  image  capture  time. 

We  have  provided  a  simple  connection  between  reflective  ghost  imaging  and  compressive 
sensing,  but  some  theoretical  issues  require  further  exploration.  If  the  target  has  a  sparse 
representation,  then  a  speckled  version  of  the  target,  should  still  have  a  sparse  representation 
(though  not  as  sparse).  Therefore,  we  should  be  able  reconstruct  the  speckled  version  of  the 
target  with  compressive  sensing.  However,  we  would  like  to  reconstruct  the  actual  target 
without  the  speckle.  We  will  look  into  ways  of  adapting  compressive  sensing  techniques  to 
do  this,  possibly  utilizing  multiple  bucket  detectors  or  multiple  wavelength  sources.  We  will 
also  explore  ways  of  measuring  the  quality  of  these  images,  and  perform  a  more  rigorous 
comparison  to  modern  pulsed  laser  radar  systems. 


78 


Appendix  A 


Detailed  SNR  Derivation 


In  this  appendix  we  will  show  in  more  detail  the  derivations  performed  to  arrive  at  the  SNR 
expressions  for  the  pseudotherm al  ghost  imager  in  Eq  (2.28)  and  the  SPDC  imager  in  Eq 
(3.18).  These  derivations  parallel  each  other,  so  we  will  not  show  as  much  detail  in  the 
derivation  of  Eq  (3.18);  the  differences  are  in  the  sources  of  noise,  the  correlation  functions, 
and  the  complexity  of  the  subsequent  integrations.  To  complete  the  SNR  derivations,  we 
need  to  calculate  the  second  moment  of  the  correlation  function  (C2(pi)). 


A.l  Pseudothermal  Ghost  Imager 

We  will  first  perform  the  derivation  for  the  pseudothermal  ghost  imager,  and  thus  start  with 
the  second  moment 

(C2(pi))  =  -L  f  dt  f  du{ii(t)i2{t)ii(u)i2(u)).  (A.l) 

1 1  J —J'i  /2  J-T,/ 2 

Using  Eq  (2.2)  to  expand  the  current  expressions,  recognizing  the  statistical  independence 
of  the  shot  noise  and  the  fields,  and  using  our  shot  noise  assumptions  from  Chapter  2  we 
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arrive  at 


J cLti  J dr2  J dr3  J dr 4  h(t  —  Ti)h(t  —  r2)h(u  —  r3)h(u  —  r4)  x 
[dV{P(p)P(p)Pi(bi)P2(p))  +  gV(Pi(r1)Pi (r3) Ai2(r2) Ai2(r4))  - 
g2r/2{P2(7'2)P2(^4)A?'i(r1)A?:1(r3)}  +  (Aii{jft)  Aii(r3)Ai2(r2)  A*2(t4))] 


Using  iterated  expectations  and  (Azm(fi)Aim(f2))  —  Q2vf>m{ti)d(t\  —  f2),  Eq  (A.l)  can  be 
written  as  the  sum  of  four  terms  as 


{CHp1))  =  T1+T2  +  T3  +  Ti 


(A. 2) 


where 


Ti  —  dr  1  /  dr2  /  dr3  /  dr4  /  dp7  /  dp"  h(t  —  Ti)h(t  —  r2)h{u  —  rA)h[u  —  r4) 


rA2 


x  c/V  A2  ( I  Pi  (Pi,  n )  I2 1 P2  (p' ,  r2)  1 2  |Pi  (pi ,  r3 )  I2 1 /A  ip?,  r4 )  | 2) , 

P  j  dr%  J dr2  J d r3  J  dp'  h(t  —  7| ;//(•/  —  r2)h{u  —  T3)h(u  —  t2) 

xgVA2(|P1(pi,ri)|2|P1(p1:r3)|2|P2(^r2)|2), 

P3  =  /  dr!  I  dr2  I  dr4  I  dp1  I  dp"  h(t  —  T\)h{t  —  r2)h(u  —  T\)h(u  —  r4) 

J  J  J  J  A2  «/  A‘2 


x  9  V  A 1  ( |  Pi  ( pi  ri )  |  ■ - 1 P2  ( pV2 )  |  ■ “  |  P2  ( p" ,  r4 )  | 2 } , 


P 


dr4  J dr2  h(t  -  rx)h{t  -  T2)h{u  -  n)h(u  -  T2)q4rf  (\El{p1.T1)\2\E2(p' ,  r2)|2). 


These  four  terms  represent  our  four  sources  of  noise:  p  conics  from  the  randomness  of  the 
fields  at  the  detectors  due  to  the  fluctuations  of  the  pseudothermal  source  and  scattering  off 
the  target:  T2  arises  from  the  randomness  in  the  field  at  the  high  spatial-resolution  detector 
beating  with  the  shot  noise  in  the  bucket  detector;  P3  comes  from  the  randomness  in  the  field 
at  the  bucket  detector  beating  with  the  shot  noise  in  the  high  spatial-resolution  detector; 
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and  finally,  i4  is  the  beat  of  the  shot  noises  from  the  two  detectors.  The  most  complicated 
of  these  terms  is  1\,  and  in  fact  the  methods  used  to  solve  the  other  terms  are  a  subset  of 
those  used  to  solve  Tf.  Therefore,  we  will  show  the  most  detail  in  the  derivation  of  7\ .  and 
show  the  parallels  for  the  other  terms. 


A. 1.1  Source  fluctuations  noise  term,  T\ 

To  evaluate  Ti  we  want  our  averaging  of  the  fields  to  be  a  jointly  Gaussian  eighth-order 
moment,  so  that  we  can  use  Gaussian-moment  factoring.  However,  the  measured  fields 
E\{p\,t)  and  E2(p.t )  are  not  jointly  Gaussian,  due  to  the  reflection  off  of  the  target.  Since 
we  are  using  a  semiclassical  treatment,  we  directly  measure  the  fields  impinging  the  detectors, 
so  E\(p,  t)  =  E'R{pi,.t)  and  E2(p,t)  —  E'T(p.  t),  the  latter  of  which  we  can  back  propagate 
to  E's(p'.t)  by  Eq  (2.3). 

With  that,  we  find  the  field  average  from  Tj  becomes 


( |  El  ( p  i .  -i )  | 1 2 1 E  2  ( P  '  •  r2 )  | : 2 1  Ei  ( p ! ,  r3 )  | : 2 1 E2  ( p"  .  t4  )  | 1 2 ) 


x  (E%(pu  TjE'xipLTjE'sipz,  t2)E's(P3,  t2)E»{Pi,  T3)E'R(Pl.TS)E%(Pi'  n)E's(ps.  r4)> 
x  {T*(p2)l1(P3)l1*(P4)T(p5))e”^(i+|p'_P2llfe^L(L+|p'_p;i,2)W^(z'+|p"_p4|2)e^(i+|p"_p5|2). 


Wo  now  have  a  fourth-order  moment  in  the  target  surface,  and  an  eighth-order  moment  in 
the  fields  to  evaluate.  Both  terms  are  higher-order  moments  of  zero-mean  Gaussian  random 
processes,  so  we  can  apply  the  Gaussian-moment  factoring  theorem  to  express  them  in  terms 
of  second  order  moments.  Since  we  are  using  pseudothermal  light,  only  the  phase-insensitive 
moments  matter;  the  phase-sensitive  cross  correlation  between  E's{p.  t)  and  ER(p ,  t),  as  well 
as  the  phase-sensitive  autocorrelation  for  each  field,  are  zero.  The  same  is  true  for  the 
surface,  which  only  has  a  phase-insensitive  autocorrelation.  Thus, 

(r (p2)r(p3)r(p4)r(p5))  =  s(P2  -  P3)s(Pa  -  p5)a^T(p2)T(p4) 

+  h(p2  -  P,5)<5(P4  -  P3)AoT(p2)T(p4).  (A. 3) 
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giving  us 

(|B1(p1,ri)|2|£,(^',r2)|2|£1(p1,r3)|2[E2(p",T4)|2)  =  -L  Jdp2J ipiT{p2)T{pi) 

x[{^(p1,T1)E^(p1,T1)£J(P2-T,)Bi(p2:r2)£g(p1,r3)E^(p1.r3)EJ(P4,T4)£s(p4,r4)> 

+  (E'-(Pl.  T2)E'S(pi,T2)E»(p1]T3)E'R(p1,  T3)E's(Pi,  Tt)E's(p2,  T,)) 

X  (A.4) 

We  now  apply  the  Gaussian- moment  factoring  theorem  to  the  eighth-order  field  moments, 
setting  all  phase-sensitive  moments  to  zero  and  applying  our  DC-block  filter,  to  arrive  at 
a  complicated  expression  that  is  a  function  of  the  propagated  Gaussian-Schell  correlation 
function  in  Eq  (2.10).  In  this  model  the  time  and  space  arguments  arc  separable,  so  we  can 
perform  the  time  integrations  (over  u,  t.  T\.  r2.  r3.  74)  and  the  spatial  integrations  (over  p2> 
Pa  p'  p")  separately. 

Time  integrations  Each  time  integral  is  of  some  function  that  is  a  product  of  correlations 
functions  and  detector  responses,  all  of  which  are  of  a  Gaussian  form.  This  means  that  the 
integrations,  while  complicated,  are  tractable  for  a  computer  to  evaluate.  The  end  result  is 
a  function  of  exponentials  and  error  functions,  but  using  our  assumption  of  Tj/Tq  3>  1  we 
can  say  that  erf[7j/7o]  ~  1  and  exp [— (T//T0)2]  «  0.  giving  us  fairly  simple  results. 

Spatial  integrations  The  spatial  integrals  arc  more  complex,  as  the  target  surface  is 
involved,  and  not  all  of  the  propagation  terms  canceled.  However,  the  the  coorelation  func¬ 
tions  are  Gaussian,  so  with  a  few  basic  manipulations  we  can  evaluate  these  integrals.  First, 
we  apply  our  previously  stated  assumptions  that  the  the  coherence  length  at  the  target  is 
small  enough  to  resolve  all  features,  and  the  intensity  radius  is  large  enough  that  the  average 
on-target  illumination  pattern  is  approximately  constant  so  that 

I dP2  T (pTh  p'  «T(pi)  J dp2e^p'-p^^.  (A.5) 
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Second,  we  define  the  Fourier  transform  of  the  target  as  T(£)  so  that 

J dp,T(p2)e-‘J?'»"’1-^  =  T  0V  - />"))  ■  (A.6) 

Finally,  we  need  to  handle  the  integrations  at  the  bucket  detector.  For  the  first  eighth-order 
term  in  Eq  (A. 4),  the  propagation  terms  cancel,  so  we  can  simply  evaluate  J^dp'  JA  dp"  = 
A%.  However,  in  the  second  term,  the  propagation  terms  do  not  all  cancel,  and  leave  complex 
exponentials  that  will  effect  a  Fourier  transform  to  the  difference  coordinates  ^(p '  —  p"). 
This  will  leave  us  needing  to  evaluate  three  different  functions  of  these  difference  coordinates. 


First,  we  have 


r  r  P2,ki.p'-P"\2  l°- 

/  d,p  /  dp"  e  2l?  =  A2^rp 

Ja>.  Jai  PlH 


I  V  | 

du  0(v.  4a). 


(A. 7) 


where  we  have  employed  the  normalized  difference  coordinate  v  —  p^k^p'  —  p")/L,  and  are 
utilizing  the  two-circle  overlap  function  and  parameter  a  =  y/Ao/nal  as  defined  in  Chapter 
2.  There  is  no  closed  form  of  this  integral  that,  we  could  identify,  but  it  is  constant  for  a 
given  a,  so  we  have  defined  the  normalized  parameter  T  in  Eq  (2.26)  so  that 


2  i  2 1  /  n  1 2 

Pjr  ^qIp  p  1 


/  dp 1  dp"  e  27? -  .Uaf.n- 

L\,  Ja2  2 


(A.8) 


Second,  we  need  to  evaluate  terms  of  the  form 


ko  /  , 


T(pi  )J  dp  J  dp"  T  yj~(p'  -  p")J  e  *rlP  p  re  j- 


(A.9) 


To  do  this  we  will  employ  the  difference  coordinate  £  =  ^(p'  —  p")  and  the  inverse  Fourier 
transform  definition 


x(p)  m  F-i [X(0](p ) 
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to  arrive  at 


T(pi)  f  dp'  f 

J  A'2  J  A.2 

=  T(Pl)A T. 

Ko 

Since  Jr“1[T(^)J(p1)  =  T(pi)  and  Jr_1[e"pi^l"/2(^)](p1)  =  (1/V p|)e~lpi l2//2pt ,  and  pL  is  small 
enough  to  resolve  all  features,  we  have  T  1[T(^)](p1)  *  T  1  [e  P^,4|2/2($)](pi)  =  T(pi).  Also, 
the  two-circle  overlap  function  0(£,2r)  is  simply  the  convolution  of  two  circles  of  radius  r 
as 


dp" T  (hy  -  p")  j 


r-'WOKpi) » j-1[e-',r£|2/2(?)|(pi) ,  r 


O  {,2 


ka  A 


L  V  7T 


(0 


(Pi)- 


7 rr20(£.  2r)  —  circ 


(A. 10) 


where  circ(-)  is  the  unit-length  circle  function.  The  inverse  Fourier  transform  of  a  circle  is  a 
.line  function  as 


Trr2^  1  [O  ($  2 r)  (0]  (pi)  =  ?’4  1  J]^f>'^  ^ 


\Pi\r 


(A. 11) 


whose  width  is  ~  pL  and  thus  resolves  the  target.  This  leaves  us  with 


T(pi)  /  dp '  I  dp"  T  (  |V  -  p"))  e-%# iF-FTe^Pi.(P'-P")  =  r-(pj.4^.  (A. 12) 


Finally,  we  need  to  evaluate  JA  dp'  JA.,dp"  |T  (^(p'  —  p"))|'.  We  again  use  difference 
coordinates  and  the  two-circle  overlap  function  to  evaluate  this,  and  define  the  term  as 
in  Eq  (2.25)  to  arrive  at 


A 


r2 

—  4' 

Kq 


(A- 13) 


It  should  be  noted  that,  by  Parseval’s  theorem,  A'r  >  A'F.  with  equality  when  T(£)  is  very 
narrow,  in  which  case  A'F  ~  /d£|T(£)|2  —  Jdp|T(p)|2  =  A'r.  With  the  final  assumption 
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of  A'T/p2L  3>  30,  we  can  combining  all  of  these  terms  to  arrive  at  an  expression  for  the  field 
noise  of 


Ti  = 


AfAlqhf^pjTo 

IAT, 


An 


T 


V2npj 


+  Tj7'T'{P\)  + 

Io 


A'f  L 2  Tj  2kL2F 

V2npi  A^kl  To  A2ko Pl  Pl 


(A.  14) 


A.  1.2  Remaining  noise  terms 


Evaluating  the  next  three  terms  only  involves  a  subset  of  the  operations  needed  evaluate 
T\.  We  again  need  to  back  propagate  to  E's(p'.t)  with  Eq  (2.3)  to  arrive  at  functions  of 
higher-order  Gaussian  moments.  For  T2  and  T-3  wo  evaluate  sixth-order  Gaussian  moments, 
and  for  T4  we  have  a  fourth-order  Gaussian  moment;  all  of  these  can  be  expressed  in  terms  of 
sums  of  products  of  the  correlation  function  in  Eq  (2.10)  by  Gaussian-momcnt  factoring.  All 
of  the  temporal  integrations  occur  over  a  products  of  Gaussians,  and  can  again  be  evaluated 
by  mathematical  software,  such  as  MATLAB  or  Mathematica.  All  of  the  spatial  integrals 
can  be  evaluated  with  the  methods  used  for  1\.  With  this,  we  arrive  at  expression  for  these 
three  noise  terms  as 


To  = 

T3  = 
T4  = 


TTCr, 


7T  cy 


ira% 


2 P  \  qirjiA‘\ A22np2L 


LrTj 


T{pi 


2 P  \3  q^r]3 Ai  Ao2tt p'2L 


l  AT, 


AirAopi  ,  L 2 


2 P  V  q^rj1  AiAoy/A p2LA 


b 


AV2lATj 
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(A. 16) 


(A.17) 


Combining  all  of  these  terms,  and  substituting  a0  —  2L/k0pL.  we  now  can  write  the 
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second  moment  of  the  correlation  function  from  Eq  (A. 2)  as 
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The  mean  of  the  correlation  function  can  be  found  in  Eq  (2.22),  giving  ns  a  squared  mean  of 
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(A.  19) 


Combining  this  with  the  second  moment  in  Eq  (A.18),  and  substituting  the  dimensionless 
brightness  term  X  =  PTop2L/a2L,  we  arrive  at  to  SNR  of 
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which  simplifies  to  Eq  (2.28),  which  wo  reproduce  here: 
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A. 2  SPDC  Ghost  Imager 


The  derivation  for  the  SPDC  ghost  imager  closely  follows  that  of  the  pseudothermal  ghost 
imager,  so  we  will  provided  a  less  detailed  derivation.  The  deviations  come  in  the  generation 
of  the  noise  terms,  their  physical  interpretation,  and  the  cross  correlation  function  between 
the  signal  and  reference  fields. 

As  before,  the  SNR  derivation  requires  evaluation  of  the  second  moment  of  the  correlation 
function 

1  r1  i  ?  pi i.  - 

(C2{px))  =  — 7  /  dt  du  {ii(t)i2{t)h(u)i2(u)) .  (A. 22) 

1I  J-T,/ 2  J-Tj/2 

We  can  use  Eqs  (3.5)  and  to  write  the  currents  in  terms  of  the  detected  field  operators 
Ei(p,t)  and  E2{p.t)  as 


—  q4A 1  f  dr i  f  dr2  [ dr3  f  dr4  f  dp'  f  dp"  h(t  —  Ti)h(t  —  r2)h(u  —  r3)h(u  —  r4) 


'A, 


x  (E{  (pi,  Ti)Ei(pi,  t\)E\(p' .  to ) E2 ( p' .  t2)E\ (pi,Tz)Ei{pi,  t3)E.\(p" .  r4)E2(p",  r4)). 


We  now  use  the  commutator  relationships  in  Eqs  (3.1)  and  (3.2)  to  normally  order  the  held 
operators  as 
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x  (EKpi.r^EKp1,  t2)E\(p" ,  t4 ) E2 ( p> .  r2)£'i(pi.  T3)E2{p".  t4))S(t i  -  r3) 
x  (E\(pi,Ti)El(p' .  t2)E\(pi,  t3)E2(p",  t4))5(p'  -  p")S(r2  -  t4)5(tx  -  r3). 


Using  Eq  (3.6)  we  can  write  this  in  terms  of  the  propagated  field  operators  E's  and  E'R. 
Since  the  terms  are  normally  ordered,  and  the  vacuum  terms  are  zero  mean  and  commute 
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with  each  other  and  the  field  operators,  this  is  simply 


(E\{Pi:  ri)E1(p1,r1)$l(p'1  T2)E2{p' T2)E\{pi:  r3)Ey(py,  T3)ib|(p",  r4)T2(p",  r4)) 

=  il{E*(Pv  n)E}(p'.,  T2)Ei(pu  r3)£?(p":  n)E'R(Pl,  n)4(p'.  r2)E'R(pur3)E'T(p".,  r4)> 
x  ^(p^f^Cp',  l|E2(pi,  r3)4(p1;  Ti)EfR(pi,  t3)E't(P",  t4))5(P'  -  p")<5(r2  -  r4) 
x  ^{EniPi- Ti)E't{p t2)£?(p",  r4)^(p/,  r2)ElR(p1,Ti)ElT{p" :  r4))5(r1  -  r3) 
x  r{^(Pl•  ■n)£1TT(p/•r2)£:fi(Pl^3)£:r(p,^r4)}h(p,  -  p")d(r2  -  r4)(5(ri  -  r3). 


We  again  have  the  correlation  second  moment  as  a  sum  of  four  terms.  (C:{pi))  —  T\  + 
7 2  +  T3  4-  / ;.  where 
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We  now  have  four  terms  which  can  be  described  as  arising  from:  the  randomness  in  the 
fields  from  the  source  fluctuations  and  scattering;  the  source  randomness  in  the  reference  field 
mixing  with  the  fluctuations  from  the  measurement  of  the  signal  field;  the  source  randomness 
in  the  signal  field  mixing  with  the  fluctuations  from  the  measurement  of  the  reference  field; 
and  the  beating  of  the  randomness  from  the  two  quantum  measurements  of  the  Gaussian 
field  operators.  For  simplicity,  we  will  evaluate  these  terms  separately 


We  first  back  propagate  E'T(p,  t)  to  E's(p%  t)  by  Eq  (3.3),  making  our  noise  terms  a  func¬ 
tions  of  normally  ordered  Gaussian  moments,  to  which  we  can  apply  the  moment  factoring 
theorem  to  obtain  an  expression  which  is  a  sum  of  the  product  of  second  order  moments. 
In  the  classical  case,  the  signal  and  reference  fields  had  phase-insensitive  cross  and  auto 
correlations  given  by  Eq  (2.6).  For  the  SPDC  ghost  imager  we  have  a  phase-sensitive  cross 
correlation  given  by  Eq  (3.7)  and  a  phase-insensitive  autocorrelation  for  each  field  given  by 
Eq  (3.8).  These  propagate  into  the  far  field  as  Eq  (3.9)  and  (3.10).  After  performing  Gaus¬ 
sian  moment  factoring  our  terms  are  the  product  two  different  correlation  functions,  one  of 
which  is  complex.  However,  after  some  additional  manipulation  the  temporal  and  spatial 
terms  can  be  separated,  and  evaluated  as  they  were  for  the  classical  case.  Substituting  in 
our  dimensionless  brightness  term  X.  the  first  term  1\  evaluates  to 
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The  remaining  noise  terms  are  evaluated  in  the  same  manner,  giving  us 
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Combining  these  terms,  and  substituting  oq  =  2L/k0pL,  we  have  an  expression  for  the  second 


moment  of  the  correlation  function 
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The  mean  of  the  correlation  function  is  in  Eq  (3.16),  and  after  substituting  in  our  bright¬ 
ness  term  X.  we  have  a  squared  mean  of 
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(A. 23) 


Combining  this  with  our  second  moment,  and  substituting  I  =  (l/(2\/27rX)  +  1).  we  arrive 
at  the  SNR  expression  in  Eq  (3.18),  which  we  reproduce  here: 
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